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ABSTRACT 


In this analysis it is shown that for satellites in circular Earth 
orbits with altitudes of 500 kilometers to 1500 kilometers, and for 
satellites in elliptical orbits with an approximate 4000 kilometer 
height of perigee, a high degree least squares polynomial (e.g. , a 
9th or 10th degree in some cases) is required to smooth both range 
and range rate data for purposes of input to orbit determination 
programs. This is especially true for overhead passes. In order 
to circumvent this problem, functions of range and range rate are 
smoothed with lower degree least squares polynomials (e.g. , 3rd 
and 4th degree) and it is shown that under the above geometric con- 
straints the standard deviation of fit (i.e. , the standard deviation 
of the errors between the approximating function for range or range 
rate and the true function of range or range rate) can be reduced to 
levels commensurate with typical S-band tracking system resolu- 
tion which is 1 to 2 meters in range and 0. 005 meters /second in 
range rate for a 1 per second data rate. The computational time 
required for fitting functions of range and range rate for low de- 
gree polynomials is shorter than the computational time required 
to fit the range or range rate directly for higher degree polynomials . 

Also shown in the analysis are the effects of Gaussian random 
noise, biases, and periodic noise. This analysis includes num- 
erous examples applied to the 44 point data smoothing interval cur- 
rently used in much of the operational preprocessing at the Goddard 
Space Flight Center. 
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SUMMARY 

In orbit determination, range, range rate, and other types of 
tracking measurements are made at various tracking sites , then 
sent to a central control for preprocessing and input to an orbit 
determination program after some smoothing technique has been 
used to reduce the effects of random and periodic noise on the data. 
For the Goddard Range and Range Rate System, both data types 
(range and range rate) are smoothed by a 3rd’ degree least squares 
polynomial (in blocks of 44 data points) and the midpoint of the 
smoothed curve is selected as the best estimate of the measure- 
ment over the data arc for input to the orbit determination program. 

In this analysis it is shown that for satellites in circular Earth 
orbits with altitudes of 500 kilometers to 1500 kilometers and for 
satellites in elliptical orbits with an approximate 4000 kilometer 
height of perigee, a high degree least square polynomial (e.g. , 9th 
or 10th in some cases) is required to smooth both range and range 
rate data for purposes of input to orbit determination programs . 

This is especially true for overhead circular passes and for satel- 
lites in the perigee portion of elliptical orbits directly overhead 
the tracking station. In order to circumvent this problem, functions 
of the range and range rate measurements are | smoothed by lower 
degree least squares polynomials (e.g., 3rd and 4th degree), and 
it is shown that the standard deviation of the errors between the 
true range and range rate and the approximating functions of range 
and range rate can be reduced to levels commensurate with typical 
S-band tracking system resolution: 1 to 2 meters in range and 0. 005 
meters /second in range rate (sample error for a sample rate of one 
per second). 

Also shown in the analysis are the effects of Gaussian random 
noise, biases, and periodic noise. 
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SMOOTHING OF FUNCTIONS OF RANGE AND RANGE RATE 
MEASUREMENTS FROM EARTH ORBITING SATELLITES 


1.0 STATEMENT OF PROBLEM 

In scientific satellite orbit determination, practical experience has shown 
the advantages of editing, filtering, and compacting tracking data prior to usage 
in an orbit determination program. Since computer running time for an orbit 
determination program is generally proportional to the number of input observa- 
tions, it has been proven economical to fit, in the least squares sense, blocks 
of consecutive raw data points to low order polynomial functions in a separate 
preprocessor computer program. Each polynomial function evaluated at the 
midpoint can be selected as one input data value to the orbit determination pro- 
gram. At Goddard for example, 44 raw data points are used as one block and a 
3rd degree least squares polynomial calculated for both range and range rate data. 
The value of the polynomial at the midpoint is then selected as input to the orbit 
determination program. By this procedure a data compaction of 44:1 of input to 
the Goddard orbit determination programs is effected. In addition, obvious "wild" 
data points can be rejected in the polynomial function calculation by comparison 
of residuals between the function evaluated at a data point and that data point. 
Editing criteria, such as 3.5 times the root-mean-square error calculation from 
all residuals, can be used to reject individual data points that clearly are not 
noise-like errors. After rejection of "wild" data points, a "closer" polynomial 
function can be fit to the remaining unedited data points in the data block. 

Along with the advantages of preprocessing raw tracking data, experience 
has shown that pitfalls are also present. It is the intent of this document to 
illustrate a primary shortcoming of preprocessor programs, and to propose an 
improvement. This demonstration of difficulty in fitting data points in the least 
squares sense to low degree polynomials is contained within Appendix A in 
Figures A-l to A-120. In generating each of the figures, the following procedure 
was applied: 

1. Generate observations from a hypothetical ground station located on 
a spherical Earth and rotating within the orbit plane of a satellite in 
a two-body system. Fix the observation sample rate at one observa- 
tion per second. 

2. Generate the data such that zero time occurs as the satellite crosses 
station zenith. 

3. Calculate the coefficients of low order polynomials fitted to the obser- 
vation data in the least squares sense. 
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4. Plot the difference between the generated observation value and the 
calculated polynomial at each data point. (These are Figures A-l to 
A-120. ) 

These differences between the polynomial function and the unperturbed ob- 
servation or "true" values are errors caused by the insufficiency of low order 
polynomials to accurately fit the equations which define the measurements of 
range , range rate and angles . 

Four different situations were considered in Figures A-l through A-120: 

1. Figures^A-1 to A-30 reflect the polynomial fitting error for range, 
range rate, and elevation data for a circular orbit, altitude of 1000 
kilometers, where the satellite passes over the station at time equals 
zero seconds, and data points are taken at time = l, 2,3. . . , 44 seconds. 
Each data type is fitted to a one through tenth degree polynomial. In 
each figure 'a* is the nominal Earth radius plus altitude and eccentricity 
’e’ is zero. 

2. Figures A-31 to A-60 show the polynomial fitting error for range, 
range rate, and elevation data for a circular orbit, altitude of 1000 
kilometers where the satellite passes over the station at time equals 
zero seconds, and data points are taken at time = 1, 2, . . . , 200 seconds. 
Figures A-31 to A-60 differ from Figures A-l to A-30 only in the length 
of time (200 seconds compared to 44 seconds) over which data is taken. 

3. Figures A-61 to A-90 depict the polynomial error for range, range 
rate, and elevation data for an elliptic orbit, semi-major axis of 95,000 
kilometers, eccentricity of 0.89 (nominal values for IMP-G), where 
the satellite passes over the station at perigee and at time equals zero 
seconds, and data points are taken at time = 1, 2, . . . , 44 seconds. 

Each data type is fitted to a one through tenth degree polynomial. 

4. Figures A-91 to A-120 differ from Figures A-61 to A-90 only in the 
span of time over which data is taken. Figures A-91 to A-120 cover 
200 seconds of time for the assumed elliptic orbit conditions. 

The figures in Appendix A show clearly the error in polynomial smoothing 
at the time of each data point on the polynomial and also shows that the number 
of zero crossings is equal to the degree of the polynomial plus one. Kruger in 
Reference 1 also considered this topic but the pictorial illustration of error is 
lacking. The report shows efficient methods of reducing this kind of error, and 
these methods achieve a significant reduction in computer operation time when 
compared to the standard method of polynomial smoothing. 
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2.0 SMOOTHING OF FUNCTIONS OF RANGE AND RANGE RATE BY LEAST 

SQUARES POLYNOMIALS 

It has been shown previously that high degree least squares polynomials 
(e.g. , 9th and 10th degree in some cases) must be fit to range and range rate 
measurements in order to reduce the rms (root mean square) of fit to levels 
commensurate with the resolution of a typical S-band system (1 to 2 meters in 
range and approximately 0. 005 meters/second in range rate). This is particu- 
larly true for satellites in a near -Earth circular orbit at the overhead portion 
of the orbit and for satellites in an elliptical orbit around the earth at the perigee 
portion of the orbit. During these portions of the orbit, accelerations and higher 
derivatives are significant. 

The purpose of this analysis is to show that lower degree least squares 
polynomials (e.g. , 3rd and 4th degree) can be fit to functions of range and range 
rate with the result that the standard deviation of fit (the standard deviation of 
fit is defined as the standard deviation of the errors between the approximating 
function and the true function of range or range rate) is within the resolution of 
a typical S-band system. We will first treat analytically the case of a satellite 
in a circular near-Earth orbit at the overhead portion of the orbit, and secondly, 
a satellite in an elliptical orbit about the Earth simulating the orbit of IMP-G at 
the perigee portion of the orbit. In addition, we will simulate errors in range 
and range rate (random, bias, and periodic errors) and note the effects of these 
errors upon the standard deviation of fit. 

Briefly, the technique consists in squaring the range measurements and 
fitting least squares polynomials to this function, then taking the square root of 
the fitted function at various points along the polynomial. It will be demonstrated 
that the standard deviation of fit between the square root of a 4th degree poly- 
nomial fit to range squared and the true range function is within the resolution 
of a typical S-band system. For the case of range rate each of the values of 
range as calculated by taking the square root of the 4th degree polynomial fit to 
the square of the range measurements is multiplied by the range rate value at 
the proper time. It should be noted that the range measurements and the range 
rate measurements do not occur at the same instant of time, that is, the time 
at the satellite for the range measurement is different from the time at the sat- 
ellite for the range rate measurement. The appropriate time for the range, 
however, can be determined from the range polynomial fit. Polynomials are 
then fit to the range rate times the range or rR (r refers to the range rate meas- 
urement and R refers to the range measurement). Values from the desired 
points along the fitted r R polynomial are then divided by the range value (again 
properly time tagged to coincide with the time at the satellite for range rate) to 
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obtain a new function. It will be seen that for a 3rd degree least squares poly- 
nomial fit to rR the resulting function is close to the true generated range rate 
function in the sense that the standard deviation of fit between the approximating 
function and the true range rate function is within the resolution of a typical S- 
band system. 


3.0 ANALYSIS FOR RANGE 

In this section we will consider the satellite to be in a circular orbit around 
the Earth, in the equatorial plane and rotating in the same direction as Earth ro- 
tation. The geometry is indicated in Figure 1 with an equatorial tracking station 
and observations referenced to the point of closest approach: 

Us = He h s 



In Figure l,co e is the angular velocity of the Earth, to s is the angular velocity 
of the satellite, R e is the radius of the Earth, R s is the radius to the satellite, 
h s is the height of the satellite above the Earth, and R is the range of the satel- 
lite from the ground tracking station. 

From the geometry in Figure 1 we can write an expression for R and expand 
it in an infinite series as follows: 

R = {R s 2 + R e 2 - 2R s R e cos (to s - co e )t} 1/2 (1) 
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Letting k y = (r s 2 + R e 2 ) I/2 , k 2 = 2R s R e /(R s 2 + R e 2 ). and k 3 = (co s 
we have 

R = k t (l - k 2 cos k 3 t) 1/2 
which can be expanded in an infinite series 

R = k l {1 - V 2 k 2 cos k 3 t - (1 • l)/(2 • 4) k 2 2 cos 2 k 3 t 
- (1 • 1 - 3)/(2 • 4 • 6) k 2 3 cos 3 k 3 t - . . .} 


co. 


or 


R 


■4-l 


T (m+%)k 2 : 


m 


m=0 


(2 m-1) it r(m+l) 


(- 1 ) 


n (k 3 t) n 
(2 n ) \l 


im 


or 


Writing the square of R we have 

R 2 = R s 2 + R e 2 - 2R s R e cos (oj s - u> e ) t 

R 2 = kj 2 (l - k 2 cos k 3 t) 
Expanding R 2 in an infinite series we have 


R 2 = k 3 2 1 1 - k 2 


oo 

L 

m = 0 


^ 3 


2m 


(2 m)! 


S ). 


(2) 


(3) 


(4) 


(5) 


The series representation in equation (5) requires fewer terms than (3) to 
accurately approximate the function. 

In order to show this, we can choose a value of t = 200 seconds and compute 
the error in range for a 4th degree polynomial for both R and R 2 . Letting h s = 
1000 kilometers (i.e. , R s = 7378 kilometers), R e = 6378 kilometers, the period 
of the satellite is 


P = 2tt 



= 2t r 


/ (7378) 3 km 3 

398603. 2 km 3 /sec 2 


or 


P = 6307 seconds 


( 6 ) 
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where 'a* is the semi-major axis of the orbit and M e is the gravitational constant 
for the Earth. Therefore 

2tt 

w s = “p" = 0.000996 radians /second 

co e = 27 t/ 24(60) 2 = 0. 0000727 radians/second 
k 3 = ( w s ■ u e) = 0.000923 radians/second 

The true range at t = 200 seconds is from equation (1) 

R = 1612. 698 kilometers 

Using equation (3) we have for m = 4 

Rj = 2 9 54.467 kilometers 

and the error is 


AR t - R x - R = 1341. 769 kilometers 
Using equation (5) for t = 200 seconds and m = 4 results in 

R 2 = 1612. 696 kilometers 


and the error is 


(7) 


AR 2 - R 2 - R = -0. 002 kilometers (8) 

It can be seen from the above example that for the same time t = 200 seconds 
and the same degree polynomial approximation of R 2 in (5) (here R 2 is being ap- 
proximated as a 4th degree polynomial in t) and R in (3) (here R is being approxi- 
mated as a 4th degree polynomial in cos (K 3 t)) taking the square root of the poly- 
nomial in R 2 is more accurate than using the polynomial in R. The reason for 
this is that when using series expansions the R 2 function converges much more 
rapidly than does the R function. That is , a lower degree polynomial can be fit 
to R 2 than to R to attain a given computation accuracy. 

4.0 ANALYSIS FOR RANGE RATE 

Recalling equation (4) of the previous section: 

R 2 = kj 2 (1 - k 2 cos k 3 t) (4) 

we can differentiate to obtain 

2RR = kj 2 k 2 k 3 sin k 3 t (9) 


6 



and combining with (2) 


ki 2 k 2 k 3 sin k 3 t 

R — (10) 

2kj (l-k 2 cos k 3 t)’ /2 


Expanding (10) in an infinite series we have 

R = 1 Ak 1 k 2 k 3 sink 3 t {l + v4k 2 cos k 3 t + (1* 3)/(2* 4) k 2 2 cos 2 k 3 t 
+ (1 • 3 • 5)/(2 -4-6) k 2 3 cos 3 k 3 t + . . . } 

- *1 k » k 3 £ (J ^> 

m= 1 

V r(nw)k 2 " (y (k,t)»Y * 

^ r <n+1) \(=o (2r) ! / 

The function 2RR in equation (9) can be expanded in an infinite series to 
obtain 


RH - k, J ki k, £] ^0^1 (12) 

n= 1 

where the factor of two has been absorbed into k 2 ' . 

The series representation in equation (12) is much easier to work with than 
equation (11) since it converges more quickly. 

In order to show this we can choose a value of t = 200 seconds and compute 
the error in range rate for a 3rd degree polynomial for both R and RR. Letting 
h s = 1000 kilometers (i.e. , R s = 7378 kilometers), R e = 6378 kilometers, we 
have for the true range rate R from equation (10) 

R = 4948.9772 meters/second 

Using equation (11) for t = 200 seconds and m = 3 we have for the range rate 
Rj = 1742. 0195 meters/second 
and the error in range rate is 

AR t = Rj - R = -3206.9577 meters/second (13) 
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Using equation (12) for t = 200 seconds, n = 3, and an approximate value 
for range of R 2 = 1612.696 kilometers obtained from equation (5) for m = 4, we 
have for the range rate 


R 2 = 4948. 9336 meters/second 
and the error in range rate is 

AR 2 = R 2 - R = -0. 0436 meters/second (14) 

In this case it is seen that for a given degree of polynomial it is much better 
to fit to RR and extract the range rate by dividing by R (R being determined from 
a fit to R 2 ) than to fit R directly. Again the reason for this is the RR function 
converges much more rapidly when expanded in series than does the R function. 


5.0 RESULTS FROM GENERATING RANGE DATA FOR A SATELLITE IN A 

CIRCULAR EARTH ORBIT AND FOR A SATELLITE IN AN ELLIPTICAL 

EARTH ORBIT 

In all cases to be described it is assumed that the satellite and the Earth 
tracking station are in the same plane (equatorial) and both rotating in the same di- 
rection. The circular orbit is illustrated in Figure 1 and the eccentric orbit 
in Figure 2. In the case of the circular orbit the satellite is rotating with an 
angular velocity of gj s and the Earth is rotating with angular velocity of to e . 

The data are taken at one/second. It should be noted that for a sample rate of 
one/second the resolution in range is on the order of 1 to 2 meters for a typical 
S-band tracking system. The standard deviation of fit is the standard deviation 
of the errors between the approximating function of range or range rate and the 
true function of range or range rate. 

Table 1 is a comparison between fitting to R 2 and fitting to R for a circular 
overhead pass. The satellite height is 1000km and the number of points is 44. 

No random noise, bias, or periodic noise has been added to the range values. 

For a fourth degree fit to R 2 the standard deviation of fit is 0. 11 x 10" 5 meters 
while for a fourth degree fit to R the standard deviation is 0. 09 meters. It can 
be seen that fitting to R 2 for a specific polynomial degree is much better than 
fitting to R. Effects of biases caused by the squaring process were noted to be 
insignificant compared to the standard deviation shown. This will be shown in 
succeeding tables . 

Table 2 is again a comparison between fitting to R 2 and R for a circular 
1000km altitude overhead pass but with the number of points increased to 200. 
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No bias, or random, or periodic noise has been added to the range. For a 
fourth degree fit to R 2 the standard deviation of fit is 0. 008 meters while for a 
fourth degree fit to R the standard deviation is 107 . 5 meters . It should be noted 
that for a third degree fit to R 2 the standard deviation of fit is 9. 09 meters which 
is much higher than the resolution of a typical S-band tracking system. 

Table 3 is a comparison between fitting to R 2 and fitting to R for an eccen- 
tric orbit simulating the orbit of the IMP-G satellite at the perigee portion of 
the orbit. The semi-major axis is 95000km and the eccentricity is 0.89. The 
number of points is 44. No random noise, bias or periodic noise has been added 
to the range values . For a fourth degree fit to R 2 the standard deviation of fit 
is 0. 16 x 10 -5 meters while for a fourth degree fit to R the standard deviation 
of fit is 0. 19 x 10 ” 3 meters. 

Table 4 is a comparison between fitting to R 2 and fitting to R for an eccen- 
tric orbit simulating the orbit of IMP-G at the perigee portion of the orbit. The 
semi-major axis is 95000km and the eccentricity is 0.89. In this case the 
number of points is 200. Again, no random noise, bias, or periodic noise has 
been added to the range values . For a fourth degree fit to R 2 the standard de- 
viation of fit is 0. 01 meters, while for a fourth degree fit to R the standard de- 
viation is 1.40 meters. 
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Table 1 

Range Analysis 


Comparison between fitting to R 2 and fitting to R for a circular overhead pass. 
h s = 1000km, n = 44 points. No random noise or bias, no periodic noise. 


Degree of 
Polynomial Fit 

a R 2 (meters) 

cr R (meters) 

1 

0.28518453 x 10 4 

0.27952192 x 10 4 

2 

0.20395929 x 10° 

0.27590778 x 10 2 

3 

0.24890964 x 10 _1 

0.30444532 x 10 

4 

0. 10513290 x 10" 5 

0.92926701 x 10" 1 

5 

0.84968737 x 10' 7 

0.58031350 x 10' 2 

6 

0.33051808 x 10“ 9 

0.33328270 x 10" 3 

7 

0.32909997 x 10" 9 

0.11612154 x 10- 4 

8 

0.33795794 x 10~ 9 

0.11953913 x 10~ 5 

9 

0.34314998 x 10 -9 

0.17982671 x 10" 7 

10 

0.30374976 x 10' 9 

0.41618490 x 10~ 8 


Table 2 

Range Analysis 


Comparison between fitting to R 2 and fitting to R for a circular overhead pass. 
h s = 1000km, n = 200 points. No random noise or bias, no periodic noise. 


Degree of 
Polynomial Fit 

o R 2 (meters) 

a R (meters) 

1 

0.51020112 x 10 5 

0.36713759 x 10 5 

2 

0.71120520 x 10 2 

0.48267764 x 10 4 

3 

0.90921999 x 10 1 

0.17821334 x 10 3 

4 

0.77245760 x 10" 2 

0.10753328 x 10 3 

5 

0.65134129 x 10' 3 

0.24781483 x 10 2 

6 

0.39555368 x 10 -6 

0.12653316 x 10 1 

7 

0.24905413 x 10' 7 

0.87747226 x 10° 

8 

0.13722842 x 10" 9 

0.21976958 x 10° 

9 

0.26623817 x 10" 9 

0.11882227 x 10" 1 

10 

0.41170183 x 10" 6 

0.93721581 x 10 -2 


t 


v 
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Table 3 

Range Analysis 


Comparison between fitting to R 2 and fitting to R for an eccentric orbit simulating 
the orbit of IMP-G at the perigee portion of the orbit, semi-major axis = 95000 
km, eccentricity = 0.89, n = 44 points. No random noise or bias, no periodic noise. 


Degree of Polynomial Fit 

o r2 (meters) 

a R (meters) 

1 

0.87843840 x 10 3 

0.87710562 x 10 3 

2 

0.54749199 x 10 _1 

0.71638679 x 10° 

3 

0.66625852 x 10 -2 

0.86594629 x 10" 1 

4 

0.15986661 x 10' 5 

0.19494156 x 10" 3 

5 

0.12877727 x 10 -6 

0.15462630 x 10 -4 

6 

0.10486462 x 10' 8 

0.58462630 x 10* 7 

7 

0.98764667 x lO’ 9 

0.34437905 x 10~ 8 

8 

0.10706711 x 10~ 8 

0.10317410 x 10" 8 

9 

0.56929905 x 10' 8 

0.98555667 x 10~ 8 

10 

0.33281828 x 10 -6 

0.43096297 x 10" 6 


Table 4 

Range Analysis 


Comparison between fitting to R 2 and fitting to R for an eccentric orbit simulating 
the orbit of IMP-G at the perigee portion of the orbit, semi-major axis = 95000 
km, eccentricity = 0.89, n = 200 points. No random noise or bias, no periodic noise. 


Degree of Polynomial Fit 

a R 2 (meters) 

a R (meters) 

1 

0.17756301 x 10 5 

0.17236634 x 10 5 

2 

0.22282125 x 10 2 

0.27444208 x 10 3 

3 

0.27383164 x 10 1 

0.29418330 x 10 2 

4 

0.13417197 x lO" 1 

0.13969013 x 10 1 

5 

0.10714135 x 10' 2 

0.76762590 x 10' 1 

6 

0.93704336 x 10 -5 

0.75774230 x 10’ 2 

7 

0. 54130339 x 10 -6 

0.17361715 x 10" 3 

8 

0.69497875 x 10" 8 

0.41028219 x lO* 4 

9 

0.13151758 x 10’ 8 

0.26686856 x 10 -6 

10 

0.23331270 x 10" 5 

0.26377387 x 10" 6 





























































For the next three tables (Tables 5,6, and 7) the satellite is in a circular 
overhead pass at a height of 1000km. No random noise, bias or periodic noise 
has been added to the range values. Table 5 is for 50 points, Table 6 for 100 
points, and Table 7 for 200 points. The standard deviations of fit are respec- 
tively, 0.23 x 10" 5 meters, 0. 14 x 10" 3 meters, and 0. 008 meters for the 4th 
degree fit for a l/second sample rate. It can be seen that for longer data arcs 
the standard deviation of fit is larger for a fixed degree polynomial. The mean 
value of the differences is given to show that no significant biases are introduced 
by the squaring process. 

Tables 8 and 9 are also for a circular overhead pass at a height of 500 kilo- 
meters and 1500 kilometers respectively. The number of points is 200, and no 
random noise, bias, or periodic noise have been added to the range values. The 
fit is to R 2 . Comparing the two tables it can be seen that for satellites in lower 
orbits the standard deviation of fit is greater. This is true since the orbital dy- 
namics are greater at lower heights . 

Table 10 is for a satellite in a circular overhead pass, the number of points 
is 200 and the satellite height is 1000km. In this case Gaussian noise has been 
added to the range values, but no bias or periodic noise. The standard deviation 
of the random noise is 2 meters and the fit is to R 2 . 


For a sample error of 2 meters standard deviation and a 200 uncorrelated 
point sample, it can be expected that the resulting standard deviation between 

the calculated range and the true range would be no better than - times 

J/200-k-l 

the 2 meter sample error, where k is the polynomial degree. For example the 
minimum possible 4th degree error would be 0. 14 meters, but note that in Table 
10, the minimum error is 0.16 meters. This should be expected since the poly- 
nomial error is also included in the 0.16 meters. Also note that increasing the 
polynomial degree increases the standard deviation. No significant bias effects 
arecausedby the squaring process as indicatedby the mean value of the difference 
between the calculated value and the true value in Table 10. 


Table 11 is for a satellite in a circular overhead pass of 1000km. The 
number of points is 200. Here a bias of 5 meters has been added to each of the 
range values, but no random or periodic noise. The standard deviation of fit 
for a fourth degree fit to R 2 is 0. 008 meters . It should be noted that although 
the standard deviation of fit is low, the bias of 5 meters has not been removed 
from the data. This is shown in Table 11 by the calculation of the mean value of 
the errors between the calculated ranges and the true values. 

Table 12 is for a satellite in a circular overhead pass of 1000km. The 
number of points is 200. In this case a bias of 10 meters has been added to 
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Table 5 

Range Analysis 


Circular orbit, overhead pass, n = 50 points, h = 1000km, no random noise or 
bias, no periodic noise. Fit to R 2 . 


Degree of Poly- 
nomial Fit to R 2 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters 

1 

0.56232410 x 10 1 

0.36697705 x 10 4 

2 

0.43846117 x 10‘ 4 

0.33792588 x 10° 

3 

0.52249932 x 10' 6 

0.41412746 x 10' 1 

4 

-0.64028427 x 10" 11 

0.22551938 x 10“ 5 

5 

-0.15133992 x 10' 10 

0.18329509 x 10^ 

6 

-0.22409949 x 10‘ 10 

0.89572210 x 10' 10 

7 

-0.15425030 x 10' 10 

0.88103885 x 10' 10 

8 

-0.22118911 x 10' 10 

0.82792086 x 10~ 10 

9 

-0.15716068 x 10" 10 

0.12349933 x 10 -8 

10 

-0.29103830 x 10~ u 

0.22028116 x 10~ 7 


Table 6 

Range Analysis 


Circular orbit, overhead pass, n = 100 points, h s = 1000km, no random noise 
or bias, no periodic noise. Fit to R 2 . 


Degree of Poly- 
nomial of Fit R 2 


Mean Value of Differences Between 
Calculated Values and True Values 


Standard Deviation 
of Fit in Meters 


0.43412137 x 10 2 


0.14070889 x 10 5 


0.76245816 x 10“ 2 


0.50892545 x 10 1 


0.31838125 x 10' 4 


0.63581281 x 10° 


0.41197927 x 10' 8 


0.13705393 x 10” 3 


-0.48894435 x lO" 10 


0.11360777 x 10 -4 


-0.20518200 x 10" 10 


0.17534497 x 10' 8 


-0.39435690 x 10" 10 


0.17001838 x 10~ 9 


-0.40017767 x 10 


■10 


0.14958597 x 10' 9 


-0.42928150 x 10' 10 


0.43774422 x 10' 8 


10 


-0.54569682 x 10' 10 


0.12456672 x 10 -6 
































































Table 7 

Range Analysis 


Circular orbit, overhead pass, n = 200 points, h s = 1000km, no random noise 
or bias, no periodic noise. Fit to R 2 . 


Degree of Poly- 
nomial of Fit R 2 

Mean Value at Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters 

1 

-0.57311709 x 10 3 

0.51020112 x 10 5 

2 

0.65392906 x 10° 

0.71120520 x 10 2 

3 

-0.13274929 x 10 _1 

0.90921999 x 10 1 

4 

0.12320237 x 10' 5 

0.77245760 x 10~ 2 

5 

-0.78485391 x 10' 7 

0.65134129 x 10' 3 

6 

-0.65556378 x 10' 10 

0.39555368 x 10" 6 

7 

-0.33687684 x 10~ 10 

0.24905413 x 10‘ 7 

8 

-0.65629138 x 10' 10 

0.13722842 x 10‘ 9 

9 

-0.32378011 x 10" 10 

0.26623817 x 10' 9 

10 

-0.56388672 x 10" 10 

0.41170183 x 10-* 


Table 8 

Range Analysis 

Circular orbit, overhead pass, n = 200 points, h s = 500km, no random noise or 
bias, no periodic noise. Fit to R 2 . 


Degree of Poly- 
nomial of Fit R 2 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters 

1 

-0.14017085 x 10 5 

0.11479317 x 10 6 

2 

0.22011470 x 10 1 

0.16554739 x 10 3 

3 

-0.28085236 x 10° 

0.21908235 x 10 2 

4 

-0.11024715 x 10' 3 

0.22829275 x 10' 1 

5 

0.13522350 x 10‘ 5 


6 

-0.48545189 x 10' 9 

0. 14739109 x 10’ 5 

7 

0.31286618 x 10' 11 

0.94197462 x 10’ 7 

8 

-0.23937901 x 10 -10 

0.14093511 x 10' 9 

9 

-0.3.6161509 x 10" 10 

0.45868546 x 10" 9 

10 

-0. 14697434 x 10' 10 

0.83170198 x 10- 6 
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Table 9 

Range Analysis 


Circular orbit, overhead pass, n = 200 points, h s = 1500km, no random noise 
or bias, no periodic noise. Fit to R 2 . 


Degree of Poly- 
nomial Fit to R 2 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters 

1 

J 0. 73014488 x 10 2 

0.31550632 x 10 5 

2 

0.10515184 x 10° 

0.36535050 x 10 2 

3 

-0.87464878 x 10‘ 4 

0.46033781 x 10 1 

4 

0.22908323 x 10' 6 

0.31902773 x 10' 2 

5 

-0.23760367 x 10' 8 

0.26649971 x 10" 3 

6 

-0.43073669 x 10' 10 

0.13167802 x 10 -6 

7 

-0.61700121 x 10' 10 

0.82462316 x 10~ 8 

8 

-0.53551048 x lO' 10 

0.16415799 x 10~ 9 

9 

-0.48894435 x 10" 10 

0.25140319 x 10“ 9 

10 

-0.44237822 x lO” 10 

0.74615831 x lO^ 6 


Table 10 
Range Analysis 


Circular orbit, overhead pass, n = 200 points, h s =1000 km, Gaussian noise with 
standard deviation of 2 meters added to range , no bias , no periodic noise . Fit to R 2 . 


Degree of Poly- 
nomial Fit to R 2 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters 

1 

-0.57326496 x 10 3 

0.51020087 x 10 5 

2 

0.50526911 x 10° 

0.71133472 x 10 2 

3 

-0.16213166 x 10° 

0.90966450 x 10 1 

4 

-0.14866744 x 10° 

0.15549463 x 10° 

5 

-0.14870137 x 10° 

0.27298950 x 10° 

6 

-0.14869614 x 10° 

0.26911890 x 10° 

7 

-0.14869552 x 10° 

0.27292880 x 10° 

8 

-0.14869539 x 10° 

0.27497664 x 10° 

9 

-0.14869537 x 10° 

0.27571468 x 10° 

10 

-0. 14869535 x 10° 

0.30533768 x 10° 
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Table 11 
Range Analysis 


Circular orbit, overhead pass, n = 200 points, h s = 1000km, no random noise, 
bias of 5 meters, no periodic noise. Fit to R 2 . 


Degree of Poly- 
nomial Fit to R 2 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters 

1 

-0.56810792 x 10 3 

0.51020242 x 10 5 

2 

0.56541091 x 10 1 

0.71140992 x 10 2 

3 

0.49867262 x 10 1 

0.90917396 x 10 1 

4 

0.50000013 x 10 1 

0.81804169 x 10- 2 

5 

0.49999993 x 10 1 

0. 54635935 x 10~ 3 

6 

0.50000000 x 10 1 

0.45810348 x 10‘ 5 

7 

0.50000000 x 10 1 

0.38591418 x 10' 5 

8 

0.50000000 x 10 1 

0.88206819 x 10- 6 

9 

0.50000000 x 10 1 

0. 59604645 x 10 -6 

10 

0.50000000 x 10 1 

0.93676052 x lO" 6 


Table 12 
Range Analysis 


Circular orbit, overhead pass, n = 200 points, h s = 1000km, no random noise, 
bias of 10 meters, no periodic noise. Fit to R 2 . 


Degree of Poly- 
nomial Fit to R 2 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters 

1 

-0.56309875 x 10 3 

0. 51020171 x 10 5 

2 

0.10654289 x 10 2 

0.71160893 x 10 2 

3 

0.99867274 x 10 1 

0.90912056 x 10 1 

4 

0.10000001 x 10 2 

0.86382979 x 10~ 2 

5 

0.99999999 x 10 1 

0.44143008 x 10' 3 

6 

0.10000000 x 10 2 

0.94133833 x 10' 5 

7 

0. 10000000 x 10 2 

0.77458523 x 10' 5 

8 

0.10000000 x 10 2 

0.17711708 x 10' 5 

9 

0.10000000 x 10 2 

0. 59604645 x 10’ 7 

10 

0.10000000 x 10 2 

0.41723251 x 10- 6 
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each of the range values, but no random or periodic noise. The standard devia- 
tion of fit for a fourth degree fit to R 2 is 0. 008 meters. The mean value of the 
errors is also shown in Table 12 and is equal to the 10 meters bias added to the 
range values. 

Table 13 is for a satellite in a circular overhead pass of 1000km height. 

The number of points is 200. In this case Gaussian random noise with a stand- 
ard deviation of 2 meters plus a bias of 10 meters have been added to the range 
values. There is no periodic noise. The standard deviation of fit for a fourth 
degree polynomial fit to R 2 is 0. 16 meters. Compare this table with Table 10. 

Table 14 is for a satellite in a circular overhead pass of 1000km height. 

The number of points is 200. Here Gaussian random noise with a standard de- 
viation of 2 meters plus a bias of 100 meters have been added to the range values. 
There is no periodic noise. The standard deviation of fit for a fourth degree 
polynomial fit to R 2 is 0. 15 meters. Compare this table with Table 13. 

Table 15 is for a satellite in a circular overhead pass of 1000km height. 

The number of points is 200. Gaussian random noise with a standard deviation 
of 2 meters plus a bias of -100 meters have been added to the range values. 

There is no periodic noise. The standard deviation of fit for a fourth degree 
polynomial fit to R 2 is 0. 16 meters. Compare this table with Tables 10, and 14. 

From the above comparison it can be seen that the squaring operation and 
taking the square root of values from the calculated polynomial has no signifi- 
cant effect on the standard deviation of fit. This is true because the biases and 
noise are very small compared to the range values. 

Table 16 is for a satellite in a circular overhead pass of 1000km height. 

The number of points is 200. No random noise or bias has been added. How- 
ever, periodic noise (i.e. , a sinusoidal variation) with one period over the 
data arc and an amplitude of 10 meters has been added to the range values. For 
a fourth degree polynomial fit to R 2 the standard deviation of fit is 7. 02 meters, 
or approximately the root mean square level of the sinusoidal variation. 

Table 17 is for a satellite in a circular overhead pass of 1000km height. 

The number of points is 200. No random noise or bias has been added. How- 
ever, periodic noise with one period over the data arc and an amplitude of 
100 meters has been added to the range values . For a fourth degree polynomial 
fit to R 2 the standard deviation of fit is 70. 19 meters. 

Comparing Tables 16 and 17 it can be seen that if the amplitude of the 
periodic noise is increased by an order of magnitude, the standard deviation of 
fit is also increased by an order of magnitude. 
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Table 13 
Range Analysis 


Circular orbit, overheadpass , n=200points, h s = 1000km, Gaussian random noise 
with standard deviation of 2 meters , bias of 10 meters , no periodic noise. Fit to R 2 . 


Degree of Poly- 
nomial Fit to R 2 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters 

1 

-0.56324662 x 10 3 

0,51019945 x 10 5 

2 

0.10505629 x 10 2 

0.71173274 x 10 2 

3 

0.98378706 x 10 1 

0.90955771 x 10 1 

4 

0.98513327 x 10 1 

0.15537387 x 10° 

5 

0.98512987 x 10* 

0.27295453 x 10° 

6 

0.98513039 x 10 1 

0.26911919 x 10° 

7 

0.98513045 x 10 1 

0.27292849 x 10° 

8 

0.98513046 x 10 1 

0.27497671 x 10° 

9 

0.98513046 x 10 1 

0.27571469 x 10° 

10 

0.98513047 x 10 1 

0.30533773 x 10° 


Table 14 
Range Analysis 


Circular orbit, overheadpass, n = 200 points, h s = 1000km, Gaussian random noise 
with standard deviation of 2 meters , bias of 100 meters, no periodic noise. Fit to R 2 . 


Degree of Poly- 
nomial Fit to R 2 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters 

1 

-0.47308153 x 10 3 

0.51018667 x 10 5 

2 

0.10050887 x 10 3 

0.71531482 x 10 2 

3 

0.99837891 x 10 2 

0.90859705 x 10 1 

4 

0.99851334 x 10 2 

0. 15454660 x 10° 

5 

0.99851299 x 10 2 

0.27264695 x 10° 

6 

0.99851304 x 10 2 

0.26912184 x 10° 

7 

0.99851304 x 10 2 

0.27292569 x 10° 

8 

0.99851305 x 10 2 

0.27497738 x 10° 

9 

0.99851305 x 10 2 

0.27571477 x 10° 

10 

0.99851305 x 10 2 

0.30533770 x 10° 
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Table 15 
Range Analysis 


Circular orbit, overhead pass, n = 200points, h s = 1000km, Gaussian random noise 
with standard deviation of 2 meters , bias of -100 meters, no periodic noise. FittoR 2 . 


Degree of Poly- 
nomial Fit to R 2 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters 

1 

-0.67344845 x 10 3 

0.51021508 x 10 s 

2 

-0.99498332 x 10 2 

0.70735443 x 10 2 

3 

-0.10016215 x 10 3 

0.91073307 x 10 1 

4 

-0.10014867 x 10 3 

0.15701463 x 10° 

5 

-0.10014870 x 10 3 

0.27334791 x 10° 

6 

-0.10014870 x 10 3 

0.26911599 x 10° 

7 

-0.10014870 x 10 3 

0.27293193 x 10° 

8 

-0.10014870 x 10 3 

0.27497589 x 10° 

9 

-0.10014870 x 10 3 

0.27571460 x 10° 

10 

-0.10014870 x 10 3 

0.30533752 x 10° 


Table 16 
Range Analysis 


Circular orbit, overhead pass, n = 200 points, h s = 1000km, no random noise 
or bias, periodic noise - one period over the data arc with amplitude of 10 
meters. Fit to R 2 . 


Degree of Poly- 
nomial Fit to R 2 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters 

1 

-0.57306064 x 10 3 

0.51018713 x 10 5 

2 

0.61348907 x 10° 

0.71388901 x 10 2 

3 

-0.12593838 x 10' 1 

0.10576680 x 10 2 

4 

-0.84866522 x 10 -4 

0.70193983 x 10 1 

5 

0.12394707 x 10' 4 

0.70555806 x 10 1 

6 

0.57059493 x 10 -6 

0.70532513 x 10 1 

7 

0.22727618 x 10" 7 

0.70533856 x 10 1 

8 

0.97767042 x 10' 9 

0.70533680 x 10 1 

9 

-0.25480404 x 10’ 9 

0.70533680 x 10 1 

10 

-0.15388650 x 10' 9 

0.70533680 x 10 1 
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Table 17 
Range Analysis 


Circular orbit, overhead pass, n = 200 points, h s = 1000km, no random noise 
or bias, periodic noise - one period over the data arc with amplitude of 
100 meters. Fit to R 2 . 


Degree of Poly- 

• 1 -n.*! JL ~ t>2 

nomiai r it tu n 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters 

i 

-0.57255206 x 10 3 

0.51004346 x 10 5 

2 

0.25027554 x 10° 

0.90946409 x 10 2 

3 

-0.65065861 x 10~ 2 

0.69946409 x 10 2 

4 

-0.84345918 x 10~ 3 

0.70192585 x 10 2 

5 

0.12241143 x 10' 3 

0.70555939 x 10 2 

6 

0. 57291405 x 10' 5 

0.70532511 x 10 2 

7 

0.21139887 x 10 -6 

0.70533855 x 10 2 

8 

0.10894728 x 10" 7 

0.70533680 x 10 2 

9 

-0.12016972 x 10- 8 

0.70533680 x 10 2 

10 

-0.15126716 x 10' 9 

0.70533680 x 10 2 
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Table 18 is for a satellite in a circular overhead pass of 1000 km height. 

The number of points is 200. No random noise or bias has been added. How- 
ever, periodic noise with only one half period over the data arc and an ampli- 
tude of 10 meters has been added to the range values. For a fourth degree poly- 
nomial fit to R 2 the standard deviation of fit is 3. 10 meters, but a bias (mean 
value) of 6.33 meters appears in the calculation because the periodic noise is 
only one half of a period. In the symmetrical cases, (Tables 16 and 17), the 
mean value is very small. In Table 18 the comparable fit to range alone (no 
squaring) is also shown. It can be seen that the squaring process has not changed 
the calculated bias that results from the periodic noise. 

Table 19 is for a satellite in a circular overhead pass of 1000km height. 

The number of points is 200. No random noise or bias has been added. How- 
ever, periodic noise with one half period over the data arc and an amplitude of 
100 meters has been added to the range values. For a fourth degree polynomial 
fit to R 2 the standard deviation of fit is 31. 04 meters, and a mean value of 63.34 
meters. This table should be compared with Table 18. 

Table 20 is for a satellite in a circular overhead pass of 1000km height. 

The number of points is 200. No random noise or bias has been added. How- 
ever, periodic noise with two periods over the data arc and an amplitude of 10 
meters has been added to the range values . For a fourth degree polynomial fit 
to R 2 the standard deviation of fit is 3.73 meters. 

Table 21 is for a satellite in a circular overhead pass of 1000km height. 

The number of points is 200. No random noise or bias has been added. Peri- 
odic noise with 4 periods over the data arc and an amplitude of 10 meters has 
been added to the range values . For a fourth degree polynomial fit to R 2 the 
standard deviation of fit is 2.30 meters. 

Table 22 is for a satellite in a circular overhead pass of 1000km height. 

The number of points is 200. No random noise or bias has been added. Peri- 
odic noise with 10 periods over the data arc and an amplitude of 10 meters has 
been added to the range values. For a fourth degree polynomial fit to R 2 the 
standard deviation of fit is 0. 98 meters. 

Table 23 is for a satellite in a circular overhead pass of 1000km height. 

The number of points is 200. No random noise or bias has been added. Peri- 
odic noise with 20 periods over the data arc and an amplitude of 10 meters has 
been added to the range values . For a fourth degree polynomial fit to R 2 the 
standard deviation of fit is 0.48 meters. 

Table 24 is for a satellite in a circular overhead pass of 1000km height. 

The number of points is 200. No random noise or bias has been added. Periodic 
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Table 18 
Range Analysis 

Circular orbit, overhead pass, n = 200 points, h s = 1000km, no random noise or 
bias , periodic noise - one half period over the data arc with amplitude of 10 meters . 


Fit to R 2 


Degree of Poly- 
nomial Fit 



Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters 

-0.56679742 x 10 3 

0.51019905 x 10 5 

0.69920845 x 10 1 

0.72002611 x 10 2 

0.63211163 x 10 1 

0.97093607 x 10 1 

0.63342420 x 10 1 

0.31036987 x 10 1 

0.63342351 x 10 1 

0.31028176 x 10 1 

0.63342352 x 10 1 

0.31028168 x 10 1 

0.63342352 x 10 1 

0.31028156 x 10 1 

0.63342352 x 10 1 

0.31028156 x 10 1 

0.63342352 x 10 1 

0.31028156 x 10 1 

0.63342352 x 10 1 

0.31028156 x 10 1 


Fit to R 


Degree of Poly- 
nomial Fit 



Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters 

0.63342352 x 10 1 

0.36713658 x 10 s 

0.63342352 x 10 1 

0.48267502 x 10 4 

0.63342352 x 10 1 

0.17823768 x 10 3 

0.63342352 x 10 1 

0.10757746 x 10 3 

0.63342352 x 10 1 

0.24974704 x 10 2 

0.63342352 x 10 1 

0.33508881 x 10 1 

0.63342352 x 10 1 

0.32245014 x 10 1 

0.63342352 x 10 1 

0.31105887 x 10 1 

0.63342352 x IO 1 

0.31028384 x 10 1 

0.63342352 x 10 1 

0.31028298 x 10 1 
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Table 19 
Range Analysis 


Circular orbit, overhead pass, n =200points, h s =1000 km, no random noise or bias, 
periodic noise - 1/2 period over the data arc with amplitude of 100 meters . Fit to R 2 . 


Degree of Poly- 
nomial Fit to R 2 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters 

1 

-0.50992007 x 10 3 

0.51018042 x 10 5 

2 

0.64035451 x 10 2 

0.85021258 x 10 2 

3 

0.63330638 x 10 2 

0.32596197 x 10 2 

4 

0.63342409 x 10 2 

0.31035110 X 10 2 

5 

0.63342351 x 10 2 


6 

0.63342352 x 10 2 

0.31028168 x 10 2 

7 

0.63342352 x 10 2 

0.31028156 x 10 2 

8 

0.63342352 x 10 2 

0.31028156 x 10 2 

9 

0.63342352 x 10 2 

0.31028156 x 10 2 

10 

0.63342352 x 10 2 

0.31028156 x 10 2 


Table 20 
Range Analysis 


Circular orbit, overhead pass, n = 200points, h s =1000km, no random noise or bias, 
periodic noise - 2 periods over the data arc with amplitude of 10 meters . Fit to R 2 . 


Degree of Poly- 
nomial Fit to R 2 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters 

1 

-0.57304880 x 10 3 

0.51019275 x 10 5 

2 

0.67885475 x 10° 

0.71363720 x 10 2 

3 

-0.11782631 x 10' 1 

0. 10294925 x 10 2 

4 

0.87714913 x 10' 3 

0.37264259 x 10 1 

5 

-0.94999505 x 10" 4 

0.67701815 x 10 1 

6 

-0.38696852 x 10" 4 

0.66901152 x 10 1 

7 

-0.13922318 x 10‘ 5 

0.70589196 x 10 1 

8 

-0.54686156 x 10" 6 

0.70419848 x 10 1 

9 

0.54413977 x 10‘ 8 

0.70539728 x 10 1 

10 

-0.10403892 x 10 -8 

0.70532474 x 10 1 


Q. ^ 
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Table 21 
Range Analysis 


Circular orbit, overhead pass, n = 200 points, h s =1000 km, no random noise or bias, 
periodic noise - 4 periods over the data arc with amplitude of 10 meters . Fit to R 2 . 


Degree of Poly- 
nomial Fit to R 2 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters 

1 

-0.57307828 x 10 3 

0.51019689 x 10 s 

2 

0.67153820 x 10° 

0.71245103 x 10 2 

3 

-0.13908271 x 10' 1 

0.98094827 x 10 1 

4 

-0.10454206 x 10‘ 3 

0.22957449 x 10 1 

5 

0.37628817 x 10' 4 

0.24174527 x 10 1 

6 

0.61536298 x 10" 4 

0.24396834 x 10 1 

7 

0.11472527 x 10~ 4 

0.34990706 x 10 1 

8 

0.97367085 x 10~ 5 

0.34483242 x 10 1 

9 

0.11239810 x 10- 4 

0.47282819 x 10 1 

10 

0.11098580 x 10" 4 

0.47616309 x 10 1 


Table 22 
Range Analysis 

Circular orbit, overhead pass, n = 200 points, h s = 1000 km, no random noise or bias, 
periodic noise - 10 periods over the data arc with amplitude of 10 meters . Fit to R 2 . 


Degree of Poly- 
nomial Fit to R 2 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters 

1 

-0.57310116 x 10 3 

0.51019944 x 10 5 

2 

0.66149021 x 10° 

0.71168460 x 10 2 

3 

-0.13717537 x 10" 1 

0.93396096 x 10 1 

4 

-0.12423038 x 10" 3 

0. 97543444 x 10° 

5 

0.30591526 x 10’ 4 

0.13130648 x 10 1 

6 

0.71848698 x 10' 5 

0.12966610 x 10 1 

7 

0.19724603 x 10' 4 

0.14410911 x 10 1 

8 

0.19955211 x 10* 4 

0.14350532 x 10 1 

9 

0.19999883 x 10‘ 4 

0.14357364 x 10 1 

10 

0.19961966 x 10' 4 

0.14401765 x 10 1 
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Table 23 
Range Analysis 


Circular orbit, overhead pass, n=200 points, h s = 1000km, no random noise or bias, 
periodic noise - 20 periods over the data arc with amplitude of 10 meters . Fit to R 2 . 


Degree of Poly- 
nomial Fit to R 2 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters 

1 

-0.57310928 x 10 3 

0.51020030 x 10 5 

2 

0.65766279 x 10° 

0.71143456 x 10 2 

3 

-0.13494873 x 10’ 1 

0.92017682 x 10 1 

4 

-0.56432316 x 10 -4 

0.47985581 x 10° 

5 

0.25893803 x 10 -4 

0.67310787 x 10° 

6 

0.11611734 x lO" 4 

0.66390283 x 10° 

7 

0.20885959 x 10" 4 

0.81419273 x 10° 

8 

0.21186552 x 10- 4 

0.80656632 x 10° 

9 

0.20405610 x 10^ 

0.89263636 x 10° 

10 

0.20430135 x 10" 4 

0.88827017 x 10° 


Table 24 
Range Analysis 

Circular orbit, overheadpass, n=200 points, h s =1000 km, no random noise or bias, 
periodic noise - 100 periods over the data arc with amplitude of 10 meters . Fit to R 2 . 


Degree of Poly- 
nomial Fit to R 2 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters 

1 

-0.57311709 x 10 3 

0.51020113 x 10 5 

2 

0.65394037 x 10° 

0.71120465 x 10 2 

3 

-0.13253378 x 10 _1 

0.90919487 x 10 1 

4 

0.22352436 x 10" 4 

0.78996626 x 10’ 2 

5 

0.20819964 x 10" 4 

0. 19801803 x lO* 2 

6 

0.20937789 x 10" 4 

0.17386356 x 10' 2 

7 

0.20910955 x 10” 4 

0.22044689 x 10" 2 

8 

0.20910024 x 10 -4 

0.21811915 x 10' 2 

9 

0.20912562 x 10" 4 

0.25762194 x 10" 2 

IQ 

0.20912495 x 10 -4 

0.25558498 x 10" 2 
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noise with 100 periods over the data are and an amplitude of 10 meters has been 
added to the range values. For a fourth degree polynomial fit to R 2 the standard 
deviation of fit is 0.008 meters. 

Summing up the above noise analysis , it is seen that purely Gaussian random 
noise added to the range values is "smoothed out" in the data processing, whereas 
biases are not. In the case of periodic noise, for a fixed data rate the higher the 
amplitude (for a fixed number of periods in the data arc) , the higher the standard 
deviation of fit. Also, for a fixed data rate the greater the number of whole pe- 
riods in the data arc (for a fixed noise amplitude), the less the standard deviation 
of fit and, in this latter case the statistics of the periodic noise approaches that 
of purely random noise. (Tables 20 through 24.) 

S-band system resolution per sample is on the order of 1 to 2 meters in 
range. For a typical circular orbit of 1000km where the satellite is passing 
directly overhead and the data arc is 200 seconds at a 1/second data rate, S- 
band system resolution is obtained with a 4th degree polynomial fit to range 
squared; whereas the same resolution is obtained with a 9th degree polynomial 
fit to the range (see Table 2). For a satellite in a circular orbit of 500km pass- 
ing directly overhead the tracking station and for a data arc of 200 seconds at 
a 1/second data rate, S-band resolution is obtained with a fourth degree fit to 
range squared (see Table 8). 


6.0 RESULTS FROM GENERATING RANGE RATE DATA FOR A SATELLITE 

IN CIRCULAR EARTH ORBIT AND FOR A SATELLITE IN AN ELLIP- 
TICAL EARTH ORBIT 

In all cases to be described it is assumed that the satellite and the Earth 
tracking station are in the same plane (equatorial) and both are rotating counter- 
clockwise. The circular orbit is illustrated in Figure 1 and the eccentric orbit 
in Figure 2. In the case of the circular orbit the satellite is rotating with an 
angular velocity of cu s and the Earth is rotating with an angular velocity of to e • 
The data are taken at l/second. It should be noted that the per sample resolu- 
tion of range rate measurements is on the order of 0. 005 meters/second for a 
typical S-band system. 

Table 25 is a comparison between fitting to rR where R is calculated from 
a 4th degree fit to R 2 , and fitting to r for a circular orbit of 1000 kilometers, 
an overhead pass with the number of points equal to 44. No random noise, bias, 
or periodic noise has been included. S-band system resolution is achieved with 
a third degree polynomial fit to rR (0. 5 x 10- 6 meters /second) whereas to ac- 
hieve the same resolution or better it takes an eighth degree polynomial fit to r . 

Table 26 is a comparison between fitting to rR and fitting to r for a circular 
orbit of 1000 kilometers , an overhead pass with the number of points equal to 
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Table 25 

Range Rate Analysis 


Comparison between fitting to rR and fitting to r for a circular orbit of 1000 kilo- 
meters, overhead pass, number of points =44, no random noise, no bias, and no 

■y 

periodic noise. Fourth order fit to R . 


Degree of Polynomial Fit 

<jf R (meters/sec) 

of (meters /sec) 

1 

0.55251670 x 10' 1 

0.74415775 x 10* 

2 

0.90169777 x 10' 2 

0.. 11006757 x 10 1 

3 

0.47860247 x 10" 6 

0.42157274 x 10’ 1 

4 

0.45710565 x 10' 7 

0.31784487 x 10' 2 

5 

0.50297534 x 10‘ 9 

0.21341051 x 10- 3 

6 

0.22455044 x lO' 10 

0.85684969 x 10' 5 

7 

0.27358737 x lO' 1 * 1 

0.99355434 x 10" 6 

8 

0. 11666926 x 10' 12 ■ 

0.16837249 x 10" 7 

9 

0.22492701 x lO' 11 

0.42796710 x lO' 8 

10 

0.15643117 x 10' 9 

0.13646136 x 10' 9 


Note: r refers to the range rate measurement and R refers to the range measurement. 

Table 26 

Range Rate Analysis 


Comparison between fitting to rR and fitting to r for a circular orbit of 1000 kilo- 
meters, overhead pass, number of points = 200. No random noise, no bias, and no 
periodic noise. Fourth order fit to R 2 . 


Degree of Polynomial Fit 

oj R (meters/sec) 

of (meters /sec) 

1 

0.41853490 x 10 1 

0.27792616 x 10 3 

2 

0.72166192 x 10° 

0.14122315 x 10 2 

3 

0.76412732 x 10' 3 

0.10873037 x 10 2 

4 

0.62547192 x 10 -4 

0.28487932 x 10 1 

5 

0.80294421 x 10' 5 

0.17028630 x 10° 

6 

0.21490185 x 10- 5 

0.14174102 x 10° 

7 

0.14087773 x 10* 6 

0.37740126 x 10' 1 

8 

0.14233777 x 10* 8 

0.22903342 x 10‘ 2 

9 

0.36281798 x 10’ 7 

0.20753311 x 10‘ 2 

10 

0.36957427 x 10‘ 8 

0.54031244 x 10‘ 3 


Note: r refers to the range rate measurement and R refers to the range measurement. 
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200. No random noise, bias, or periodic noise has been included. S-band sys- 
tem resolution is achieved again with a third degree polynomial fit to rR whereas 
to achieve the same resolution or better it takes a tenth order polynomial fit to r. 


Table 27 is a comparison between fitting to rR and fitting to r for an eccentric 
orbit simulating the orbit of IMP-G at the perigee portion of the orbit, semi- 
major axis of 95000 kilometers, eccentricity of 0.89, number of points equal to 
44 with a fourth degree polynomial fit to R 2 . S-band system resolution is ac- 
hieved with a second degree polynomial fit to rR (0. 0024 meters/second) whereas 
the same or better resolution is achieved with a third degree polynomial fit to r. 
No random noise, bias, or periodic noise are present. 


Table 28 is a comparison between fitting to rR and fitting to r for an eccen- 
tric orbit simulating the orbit of IMP-G at the perigee portion of the orbit, semi- 
major axis of 95000 kilometers, eccentricity of 0. 89, number of points equal to 
200 with a fourth degree polynomial fit to R 2 . Here S-band system resolution is 
achieved with a third degree fit to rR (0. 0013 meters /second) whereas the same 
or better resolution is achieved with a fifth degree fit to r (0. 0011 meters/ 
second). No random noise, bias, or periodic noise are present. 

Table 29 is for a satellite in an eccentric orbit simulating the orbit of IMP-G 
at the perigee portion of the orbit, semi-major axis of 95000 kilometers, eccen- 
tricity of 0. 89, number of points equal to 44 with a third degree fit to R 2 . S- 
band system resolution is achieved with a second degree polynomial fit to rR 
(0.0024 meters /second). No random noise, bias, or periodic noise are present. 
Compare this with Table 27 (4th degree fit to R 2 ). The mean value of the differ- 
ences is shown to illustrate that a bias is not introduced by the method of fitting. 


Table 30 is for a satellite in an eccentric orbit simulating the orbit of IMP-G 
at the perigee portion of the orbit, semi -major axis of 95000 kilometers, eccen- 
tricity of 0. 89, number of points equal to 200 with a third degree fit to R 2 . No 
random noise, bias, or periodic noise is present. The standard deviation of fit 
for a third degree polynomial fit to rR is 0.0021 meters /second, within the re- 
solution of a typical S-band system . Compare this with Table 2 8 (4th degree fit to R 2 ) . 


Table 31 is for a satellite in an eccentric orbit simulating the orbit of IMP-Gr 
at the perigee portion of the orbit, semi-major axis of 95000 kilometers, ec- 
centricity of 0. 89, number of points equal to 200. Guassian random noise with a 
standard deviation of 0. 005 meters /second has been added to the range rate val- 
ues. No bias or periodic noise has been added and a fourth degree polynomial 
has been fit to R 2 . For a third degree fit to rR the standard deviation of fit is 
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Table 27 

Range Rate Analysis 

Comparison between fitting to rR and fitting to r for an eccentric orbit simulating the 
orbit of IMP-G at the perigee portion of the orbit, semi-major axis = 95000km, eccen- 
tricity = 0. 89, number of points = 44. Fourth order fit to R 2 . 


Degree of Polynomial Fit 

aj R (meters/sec) 

of (meters/sec) 

1 

0.14843337 x 10' 1 

0.19416260 x 10° 

2 

0.24136823 x 10' 2 

0.31365064 x 10' 1 

3 

0.72864384 x 10' 6 

0.88834831 x HT 4 

4 

0.70570925 x 10' 7 

0.84818951 x 10* 5 

5 

0.95270103 x 10’ 10 

0.37503272 x 10' 7 

6 

0.64165907 x 10' 11 

0.25111475 x 10' 8 

7 

0.21837252 x 10‘ 12 

0.22702691 x 10' 10 

8 

0.20768250 x 10' 12 

0.83556286 x 10~ n 

9 

0.37381279 x 10* 12 

0.70736022 x 10' 11 

10 

0.33854018 x lO’ 10 

0.43032299 x 10' 10 


Note: r refers to the range rate measurement and R refers to the range measurement. 

Table 28 

Range Rate Analysis 


Comparison between fitting to rR and fitting to r for an eccentric orbit simulating the 
orbit of EMP-G at the perigee portion of the orbit, semi-major axis = 95000km, eccen- 
tricity = 0.89, number of points = 200. Fourth order fit to R 2 . 


Degree of Polynomial Fit 

Of R (meters /sec) 

a- (meters /sec) 

1 

0.13235152 x 10 1 

0.16209897 x 10 2 

2 

0.21728503 x 10° 

0.23296691 x 10 1 

3 

0.13335383 x 10' 2 

0. 13825187 x 10° 

4 

0.12431866 x 10' 3 

0.91762163 x 10' 2 

5 

0.34060858 x 10' 5 

0.10524417 x 10* 2 

6 

0.14739791 x 10- 6 

0.27904480 x 10"* 

7 

0.20301123 x 10' 7 

0.73389407 x 10-5 

8 

0.14609493 x 10' 9 


9 

0.14172266 x 10' 9 

0.46887330 x 10' 8 

10 

0.19497726 x 10~ 8 

0.21698210 x 10’ 8 


Note: r refers to the range rate measurement and R refers to the range measurement. 
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Table 29 

Range Rate Analysis 


Eccentric orbit simulating the orbit of IMP-G at the perigee portion of the orbit, semi- 
major axis = 95000km, eccentricity = 0.89, number of points = 44. No random noise, 
bias or periodic noise. Third order fit to R 2 . 


Degree of Poly- 
nomial Fit to rR 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters /Sec 

1 

-0.31415660 x 10~ 5 

0.14843303 x 10' 1 

2 

0.11973341 x 10‘ 8 

0.24134640 x 10' 2 

3 

-0.94828186 x lO" 13 

0.12075541 x 10' 5 

4 

-0.40391932 x 10' 13 

0.28866104 x 10" 6 

5 

-0.31939098 x 10~ 13 

0.60276757 x 10' 9 

6 

-0.27361951 x 10~ !3 

0.86537384 x 10' 10 

7 

-0.21573652 x lO’’ 3 

0.36686199 x 10’ 12 

8 

-0.25873243 x 10' 13 

0.20022048 x 10‘ 12 

9 

-0.24621719 x 10~ 13 

0.40106217 x 10' 12 

10 

-0.48557118 x lO’ 13 

0.46081116 x 10' 10 


Note: r refers to the range rate measurement and R refers to the range measurement. 


Table 30 

Range Rate Analysis 


Eccentric orbit simulating the orbit of IMP-G at the perigee portion of the orbit, semi- 
major axis = 95000km, eccentricity = 0.89, number of points = 200. No random noise, 
bias or periodic noise. Third order fit to R 2 . 


Degree of Poly- 
nomial Fit to r R 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters /Sec 

1 

-0.50857727 x 10~ 2 

0.13234910 x 10 1 

2 

0.39501185 x 10' 4 

0.21694633 x 10° 

3 

-0.34286505 x 10' 7 

0.21404468 x 10~ 2 

4 

0.74316297 x 10' 9 

0.49350656 x 10' 3 

5 

-0.13638624 x lO' 11 

0.20279421 x 10~ 4 

6 

-0.14587997 x lO’ 12 

0.24368930 x 10' 5 

7 

-0.22972513 x 10' 12 

0.21041812 x lO” 6 

8 

-0.19679591 x 10' 12 

0.11313188 x lO" 7 

9 

-0.16329160 x 10~ 12 

0.18525224 x 10‘ 8 

10 

-0.24185765 x lO’ 12 

0.21254859 x 10~ 8 


Note, r refers to the range rate measurement and R refers to the range measurement. 
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Table 31 

Range Rate Analysis 

Eccentric orbit simulating the orbit of IMP-G at the perigee portion of the orbit, semi- 
major axis = 95000km, eccentricity = 0.89, number of points = 200. Gaussian/ noise with 
standard deviation of 0.005 meters/second added to range rate measurements. No bias 
or periodic noise. Fourth degree fit to R 2 . 


Degree of Polynomial 
Fit to rR 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation of 
Fit in Meters/Second 

1 

-0.52521228 x 10~ 2 

0.13235272 x 10' 

2 

-0.12710253 x 10- 3 

0.21728637 x 10° 

3 

-0.16666557 x 10‘ 3 

0.14042741 x 10- 2 

4 

-0.16663907 x 10~ 3 

0.48046468 x lO’ 3 

5 

-0.16663930 x 10' 3 

0.46270123 x 10‘ 3 

6 

-0.16663933 x 10‘ 3 

0.60050733 x 10‘ 3 

7 

-0.16663933 x 10' 3 

0.74288198 x 10~ 3 

8 

-0.16663933 x 10~ 3 

0.76385666 x 10~ 3 

9 

-0.16663933 x 10' 3 

0.78017102 x 10' 3 

10 

-0.16663933 x 10' 3 

0.79971799 x 10‘ 3 


Note: r refers to the range rate measurement and R refers to the range measurement. 


0. 0014 meters /second. Compare Table 31 with Table 28 to see the effects of 
the Gaussian random noise. In this case the r noise is small compared to the 
value of r. 

Table 32 is for a satellite in a circular orbit of 500 kilometers , an over- 
head pass, the number of points equal to 200, and a fourth degree fit to R 2 . No 
random noise, bias, or periodic noise has been added. S-band system resolu- 
tion is achieved with a third degree polynomial fit to rR where the standard de- 
viation of fit is 0. 0023 meters/second. 

Table 33 is for a satellite in a circular orbit of 1500 kilometers, an over- 
head pass, the number of points equal to 200, and a fourth degree polynomial fit 
to R 2 . No random noise, bias, or periodic noise has been added. The standard 
deviation of fit for a third degree polynomial fit is 0.32 x 10~ 3 meters /second. 

Table 34 is for a satellite in a circular orbit of 1000 kilometers, an over- 
head pass, the number of points equal to 200, and a fourth degree polynomial fit 
to R 2 . No random noise or bias has been added. Periodic noise has been added 
with a range amplitude of 2 meters and a range rate amplitude of 1 meter/ 
second. There are two periods over the data arc. The standard deviation of 
fit for a 3rd degree polynomial fit to rR is 0. 36 meters/second which is caused 
by the r cyclic noise. 




































Table 32 

Range Rate Analysis 


Circular orbit of 500km, overhead pass, number of points = 200. No random noise 
or bias, no periodic noise. 4th degree fit to R 2 . 


Degree of Poly- 
nomial Fit torR 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters /Sec 

1 

0.47953194 x 10° 

0.96246538 x 10 1 

2 

0.24847058 x 10' 1 

0.17378154 x 10 1 

3 

0.29492329 x 10"* 

0.22904419 x 10' 2 

4 

-0.59403423 x 10 -6 

0.14613763 x 10' 3 

5 

-0.52965080 x 10" 8 

0.46734438 x 10"* 

6 

0.10620666 x 10‘ 7 

0.32032858 x 10"* 

7 

0.29504444 x 10’ 8 

0. 11446048 x 10"* 

8 

0. 17628773 x 10‘ 9 

0.19062020 x 10" 5 

9 

-0.32519232 x 10’ 11 

0.13304742 x 10' s 

10 

0.57038285 x lO' 11 

0.89235640 x 10" 6 


Note: r refers to range rate measurement and R refers to the range measurement. 


Table 33 

Range Rate Analysis 


Circular orbit of 1500km, overhead pass. Number of Points = 200. No random 
noise or bias, no periodic noise. 4th degree fit to R 2 . 


Degree of Poly- 
nomial Fit to rR 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters /Sec 

1 

-0.25329401 x 10' 1 

0.21625184 x 10 1 

2 

0. 10486964 x 10' 2 

0.36537740 x 10° 

3 

0.76490996 x 10’ 8 

0.31597987 x 10' 3 

4 

0.20410022 x 10 -8 

0.28381520 x 10"* 

5 

0.15667219 x 10 _1 ° 

0.18873398 x 10' 5 

6 

-0.57786664 x 10' 12 

0.18084428 x 10^ 

7 

-0.41014303 x lO' 12 

0.39209944 x 10’ 7 

8 

-0.42284398 x 10* 12 

0.10342901 x lO' 7 

9 

-0.27872815 x 10" 12 


10 

-0.49690030 x 10‘ 12 

0.38911794.x 10‘ 8 


Note, r refers to the range rate measurement and R refers to the range measurement. 
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Table 35 is for a satellite in a circular orbit of 1000 kilometers, an over- 
head pass, the number of points equal to 200, and a fourth degree polynomial fit 
to R 2 . No random noise or bias has been added. However, periodic noise has 
been added to the data with a range amplitude of 2 meters and a range rate am- 
plitude of 1 meter /second. There is one period over the data arc. The standard 
deviation of fit for a 3rd degree polynomial fit to rR is 0.71 meters/second. 

Table 36 is for a satellite in a circular orbit of 1000 kilometers, an over- 
head pass, the number of points equal to 200, and a fourth degree polynomial fit 
to R 2 . No bias or periodic noise has been added. However, Gaussian random 
noise has been added with a range standard deviation of 2 meters and a range 
rate standard deviation of 0.005 meters /second. For a 3rd degree polynomial 
fit to rR the standard deviation of fit is 0. 71 x 10 -3 meters/second. 

Table 37 is for a satellite in a circular orbit of 1000 kilometers, an over- 
head pass, the number of points equal to 200, and a fourth degree polynomial fit 
to R 2 . Gaussian random noise has been added with a range standard deviation of 
2 meters and a range rate standard deviation of 0.005 meters /second. A bias 
of 5 meters has been added to the range measurements and a bias of 5 meters/ 
second added to the range rate measurements . Periodic noise has also been 
added to the data with a range amplitude of 2 meters and a range rate amplitude 
of 1 meter/second. There are two periods over the data arc. For a 3rd degree 
polynomial fit to rR the standard deviation of fit is 0. 36 meters /second. Com- 
pare this with Table 34, the noise free case, and note that the mean value con- 
verges to the bias of 5 meters /second. 

Table 38 is for a satellite in a circular orbit of 1000 kilometers, an over- 
head pass, the number of points equal to 200, and a fourth degree polynomial fit 
to R 2 . Gaussian random noise with a range standard deviation of 2 meters and 
a range rate standard deviation of 0.005 meters/second has been added to the 
data. No biases have been added. However, periodic noise with a range am- 
plitude of 2 meters and a range rate amplitude of 1 meter/second has been added. 
There are two periods over the data arc. For a 3rd degree polynomial fit to 
rR the standard deviation of fit is 0. 36 meters /second. Compare with Tables 
34 and 37. 

Summing up the previous range rate analysis it is seen that for satellites 
directly overhead the tracking station in circular orbits of 500 kilometers and 
greater and for satellites that are in elliptical orbits at the perigee portion of 
the orbit (where the perigee height is only a few thousand kilometers), in the 
noise free case where the data arcs are 200 seconds or less (at a sample rate 
of one /second) typical S-band resolution (on the order of 0.005 meters /second) 
can be achieved by fitting a 3rd degree least squares polynomial to rR (and a 4th 
degree least squares- polynomial to R 2 ) (see Tables 25, 26, 27, 28, 32, and 33). 
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Table 34 

Range Rate Analysis 

Circular orbit of 1000km, overhead pass. Number of points = 200. No random noise or 
bias. Periodic noise-range amplitude of 2 meters, range rate amplitude of 1 meter/second: 
two periods over the data arc. 4th degree fit to R 2 . 


Degree of Poly- 
nomial Fit to rR 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters/Sec 

1 

-0.40882026 x 10"' 

0.42674498 x 10 1 

2 

0.91430595 x 10' 2 

0.79177834 x 10° 

3 

0.14924478 x 10~ 3 

0.36361170 x 10° 

4 

0.86128632 x 10- 4 

0.37262966 x 10° 

5 

-0.96187081 x 10' s 

0.67699432 x 10° 

6 

-0.40732068 x 10' 5 

0.66901158 x 10° 

7 

-0.12848551 x 10" 6 

0.70589165 x 10° 

8 

-0.60574703 x 10* 7 

0.70419849 x 10° 

9 

0.12498209 x 10~ 8 

0.70539724 x 10° 

10 

-0. 13428403 x 10‘ 9 

0.70532474 x 10° 


Note: f refers to the range rate measurement and R refers to the range measurement. 


Table 35 

Range Rate Analysis 


Circular orbit of 1000km, overhead pass, number of points = 200. No random noise or 
bias. Periodic noise - range amplitude of 2 meters, range rate amplitude of 1 meter/ 
second: one period over the data arc. 4th degree fit to R 2 . 


Degree of Poly- 
nomial Fit to rR 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters /Sec 

1 

-0.46776306 x 10' 1 

0.43560247 x 10 1 

2 

0.25762243 x 10' 2 

0.90776251 x 10° 

3 

0.69169003 x 10' 4 

0:70877599 x 10° 

4 

-0.86556339 x 10' 5 

0.70192154 x 10° 

5 

0. 12504086 x 10' 5 

0.70556025 x 10° 

6 

0.57350466 x 10' 7 

0.70532509 x 10° 

7 

0.22965466 x 10' 8 

0.70533855 x 10° 

8 

0.11151414 x 10' 9 

0.70533680 x 10° 

9 

-0.10936816 x 10' 10 

0.70533680 x 10° 

10 

-0.65533357 x 10' 12 

0.70533680 x 10° 


Note: f refers to the range rate measurement and R refers to the range measurement. 
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Table 36 

Range Rate Analysis 


Circular orbit of 1000km, overhead pass, number of points = 200. Gaussian random noise- 
range standard deviation of 2 meters and range rate standard deviation of 0.005 meters/ 
second. No bias or periodic noise. 4th degree fit to R 2 . 


Degree of Polynomial 
Fit to rR 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation of 
Fit in Meters /Second 

1 

-0.42710922 x 10" 1 

0.41855463 x 10' 1 

2 

0.64746148 x 10' 2 

0.72173969 x 10° 

3 

-0.16561654 x 10‘ 3 

0.70583749 x 10~ 3 

4 

-0.16668431 x 10' 3 

0.46814769 x 10~ 3 

5 

-0.16669457 x 10‘ 3 

0.46984210 x 10' 3 

6 

-0.16664764 x 10' 3 

0.57422594 x 10‘ 3 

7 

-0.16663969 x 10~ 3 

0.76158320 x 10' 3 

8 

-0.16663899 x 10' 3 

0.77240557 x 10' 3 

9 

-0.16663938 x 10‘ 3 

0.78505025 x 10' 3 

10 

-0.16663933 x 10‘ 3 

0.80063179 x 10' 3 


Note: r refers to the range rate measurement and R refers to the range measurement. 


Table 37 

Range Rate Analysis 

Circular orbit of 1000km, overhead pass, number of points = 200. Gaussian random noise: 
range standard deviation of 2 meters; range rate standard deviation of 0.005 meters /second. 
Range bias of 5 meters: range rate bias of 5 meters /second. Periodic noise: range ampli- 
tude of 2 meters; range rate amplitude of 1 meter /second. Two periods over the data arc. 
4th degree fit to R 2 . 


Degree of Polynomial 
Fit to rR 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation of 
Fit in Meters /Second 

1 

0.49608727 x 10 1 

0.41174103 x 10 1 

2 

0.50091630 x 10* 

0.81110932 x 10° 

3 

0.49999835 x 10 1 

0.36369805 x 10° 

4 

0.49999195 x 10 1 

0.37283249 x 10° 

5 

0.49998237 x 10 1 

0.67711807 x 10° 

6 

0.49998293 x 10 1 

0.66909882 x 10° 

7 

0.49998332 x 10 1 

0.70613452 x 10° 

8 

0.49998333 x 10 1 

0.70442591 x 10° 

9 

0.49998334 x 10 1 

0.70561445 x 10° 

10 

0.49998334 x 10 1 

0.70554436 x 10° 


Note: f refers to the range rate measurement and R refers to the range measurement. 
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Table 38 

Range Rate Analysis 


Circular orbit of 1000km, overhead pass, number of points = 200. Gaussian 
random noise: range amplitude of 2 meters; range rate amplitude of 0.005 
meters /second. No biases. Periodic noise: range amplitude of 2 meters; 
range rate amplitude of 1 meter/second. Two periods over the data arc. 4th 
degree fit to R 2 . 


Degree of Poly- 
nomial Fit to rR 

Mean Value of Differences Between 
Calculated Values and True Values 

Standard Deviation 
of Fit in Meters /Sec 

1 

-0.41050371 x 10' 1 

0.42676703 x 10 1 

2 

0.89757266 x 10' 2 

0.79197794 x 10° 

3 

-0.17495436 x 10 -4 

0.36379782 x 10° 

4 

-0.80544685 x 10- 4 

0.37282381 x 10° 

5 

-0 . 17631361 x 10' 3 

0.67713959 x 10° 

6 

-0.17072026 x 10“ 3 

0.66909879 x 10° 

7 

-0.16676820 x 10' 3 

0.70613482 x 10° 

8 

-0.16669956 x 10' 3 

0.70442590 x 10° 

9 

-0.16663813 x 10' 3 

0.70561445 x 10° 

10 

-0.16663946 x 10' 3 

0.70554436 x 10° 


Note: r refers to the range rate measurement and R refers to the range measurement. 


In the case of a satellite in a 1000 kilometer circular orbit at the overhead 
portion of the orbit for a 200 second data arc where the data rate is l/second and 
there is no random noise, bias, or periodic noise on the range rate data, the 
standard deviation of fit (between the approximating function for range rate and 
the true range rate function) is 0.76 x lb' 3 meters /second for a third degree 
least squares polynomial fit to the range times the range rate. For a third de- 
gree polynomial fit to the range rate the standard deviation of fit (between the 
fitted range rate polynomial and the true range rate function) is 10. 9 meters/ 
$econd. In order to achieve a standard deviation of 0. 76 x 10' 3 meters /second 
for a least squares polynomial fit to range rate directly, a 10th degree polynomial 
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is required (see Table 26). The advantage of a lower degree polynomial fit to 
the range rate times the range occurs because the computation time for fitting 
a 3rd degree polynomial to this mixed data type is much less than the computa- 
tion time required to fit a 10th degree polynomial to range rate directly. This 
will be shown in Section 8. 

When purely random Gaussian noise with a standard deviation of 0. 005 
meters/second is added to the true range rate, the noise is "smoothed out" in 
the fitting process (compare Table 28 with 31 and Table 26 with 36). 

In the case of periodic noise, the greater the number of whole periods in 
the data arc (for the same noise amplitude), the smaller the standard deviation 
of fit (compare Tables 34 and 35). 

The effect of adding a bias to the true range rate values is not removed in 
the smoothing process, but shows up as the mean value of the differences be- 
tween the approximating function for range rate and the true function of range 
rate (see Table 37). 


7.0 COMPARISON OF RANGE RATE POLYNOMIAL ERROR FOR VARIOUS 
METHODS OF SMOOTHING — 200 SECOND DATA ARC, CIRCULAR 
ORBIT OF 1000 KILOMETER HEIGHT 

In each of the preceeding tables the effects of polynomial degree on the poly- 
nomial smoothing error for one specific case were shown. In this section a 
specific orbit is assumed: 200 second length of data with 200 data points for a 
1000 kilometer height, circular orbit, passing directly over the station at time 
equal to zero seconds; and comparisons are made for different methods of 
smoothing. The data has been simulated for a 2 dimensional situation as before. 

In Table 39, no data noise has been assumed so that the data are "perfect". 
In column 1, a polynomial is fitted to data points consisting of the true range 
(R x ) and the true range rate (R x ) multiplied together (R t Rt). The values of the 
calculated range rate (Rcalc) are obtained by dividing the polynomial value at 
t = 1, 2, 3 . . . , by the true range value. The error is formed by subtracting 
the true range rate from the calculated range rate (error = Rcalc - r t)* Biases 
and standard deviations were computed and summed in a root-sum-square man- 
ner. The standard deviation of these calculated errors are shown in column 1. 
Note that the computer truncation error introduces a minimum error point at an 
8th degree polynomial. For higher degree polynomials, an infinite digit computer 
would show smaller error. Here the error increased for a 9th degree polyno- 
mial. The specific computer used here works with 64 bits or approximately 16 
significant digit s. 
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Range rate polynomial smoothing error for various methods of fitting a 200 second arc, 1000 kilometer height, 
circular orbit, no noise on data. 





In column 2, the calculated range, Rcalc » * s f° und by squaring the true 
range and fitting the squared data points to a 6th degree polynomial in range 
squared. Polynomials are then fitted to values of the true range rate (R x ) and 
the calculated range (Rcalc) multiplied together. The error is defined as above. 
Note that the minimum error is larger than for column 1. Minimum errors 
are shown enclosed in a rectangle. 


Columns 3 and 4 are similar to column 2 except that the range values cal- 
culated (R CAL c) are formed from a 4th and 2nd degree polynomial to Rt squared. 

Column 5 shows the error from a fit to the true range rate and is the poorest 
method of fitting as stated previously , since low degree polynomials are poor 
approximations to the true range rate function. 


Table 39 is an artificial case since the true data are always perturbed by 
noise. The noise-free conditions of Table 39 were presented only for illustra- 
tion of polynomial smoothing error. Table 40 shows the more realistic com- 
parisons for data with typical noise-like errors. Column 1 is also artificial in 
the sense that it includes range noise only, that is, range rate data is assumed 
perfect. The noisy range data as squared, fit to a 4th degree polynomial, 
square root values are taken from the calculated polynomial (R CA lc ) » and the 
calculated range and the true range rate (multiplied together) are fit to various 
degree polynomials. The calculated range rate, Rcalc » results from dividing 
values from the polynomial by R CALC • A minimum occurs at a 10th degree 
polynomial. 


Column 2 uses a value of R C alc from squaring R x and fitting to a 4th degree 
polynomial. Column 2 is also unrealistic since there is no noise on the range 
data. Noise with standard deviation of 0. 005 meters per second is imposed upon 
R x , and R calc and the noisy range rate data are multiplied together, Rcalc 
(R x + N), and fit to a polynomial of degree 1 through 10. A minimum value 
of 0. 0005 meters/second occurs for a 5th degree polynomial fit. Here the 


minimum value of error should be yg' oo k 1 times va l ue °f the sample stand- 
ard deviation (0.005 meters /second) or 0.00036 meters/second for a 5th degree 
polynomial fit (k = 5). 


Columns 3, 4, and 5 have noise on both the range and the range rate data; 
the differences are in the method of calculating the range (Rcalc)* Column 3 
calculates range from a 6th degree polynomial to range squared, column 4 from 
a 4th degree polynomial, and column 5 from a 2nd degree polynomial. Although 
column 3 fits to a higher degree polynomial in (R x + N) 2 , computer truncation 
error seems to cause a larger minimum error than for the 4th degree fit to 
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(Rt+ N) 2 in column 4. In column 5, the 2nd degree fit to (Rt + N) 2 is not suf- 
ficient to calculate any useful value of range (Rcalc)* In column 5, the minimum 
error occurs for an 8th degree polynomial. 


Column 6 extracts a calculated range from a 4th degree fit to R x + N. Note 
that this is not a fit to the square of the data. Here again the calculated range, 
R C alc* is poor and minimum error occurs for a 10th degree polynomial. 

Column 7 are errors for a polynomial fit directly to noisy range rate. 
Minimum error occurs for a 10th degree polynomial fit and is approximately 
twice the value for column 4 fit to 4th degree. 


8.0 TIMING OF METHODS 


The contention has been made that performing operations on the data points 
(such as squaring the range) prior to polynomial fitting in the least squares 
senses, reduces the degree of the polynomial required, and thus conserves com- 
puter operation time. The purpose of this section is to describe some timing 
experiments made which demonstrate the validity of this claim. 


Several programs were run on the IBM 360-75 computer in order to time 
various operations. These runs were made with no other jobs in the computer 
so that complications in timing would not be encountered. In addition the ma- 
chine timing function, TIME , was called before and after each operation to be 
timed. TIME returns the number of hundredths of seconds at the time of the 
call which have occurred since midnight on that day. Since computer operations 
take relatively small time intervals compared to hundredths of a second, the 
identical operation was repeated many times between calls to TIME , in order to 
calculate an accurate estimate for a single operation. For example, the square 
root operation takes only 73. 5 /x seconds on the 360/75 computer but timing was 
performed on 20,000 of these operations. In timing the polynomial least squares 
fitting, the program used was D00107 from the GSFC computer library. In pro- 
gram D00107 the normal equations are formed, and the pivotal element method 
is used for the matrix inversion process. Table 41 illustrates the results of 
this timing. 
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Table 41 

Timing of Methods 


METHOD 1: 


Square 200 Range Measurements 

0. 0039 sec 

Generate 4th Degree Polynomial fit to 200 
Squared Measurements 

0.4140 sec 

Square Root 200 Values from Squared 
Polynomial 

0. 0160 sec 

TOTAL 

0.4339 sec 

METHOD 2: 


Generate 9th Degree Polynomial fit to 200 
Range Measurements 

1. 0796 sec 

R VTIO" METHOD 1 - o 5 
RATI °- METHOD 2 2 * * ‘ 5 



9.0 CONCLUSIONS 

In this analysis it has been shown that: 

1. High degree least squares polynomials (e.g. , 9th and 10th degree in 
some cases) are required to fit range and range rate measurements in 
order to reduce the root mean square error of fit to levels that are com- 
mensurate with the resolution of a typical S-band tracking system: 1 to 

2 meters in range and 0. 005 meters /second in range rate (sample error 
for sample rate of one per second). This is especially true for satel- 
lites that are in a near-Earth circular orbit or satellites that are in an 
elliptical orbit at the perigee portion of the orbit, where perigee height is 
only a few thousand kilometers . 

2. A method to circumvent the problem above is to least squares poly- 

nomial fit functions of range and range rate, and extract the range and 

range rate from these functions. For a given accuracy the degree of 

the least squares polynomial required to fit these functions is much 
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lower (3rd and 4th degree) than the degree of the least squares poly- 
nomial required to fit the range and range rate directly (9th and 10th 
degree in some cases). 

3. By fitting to the square of range, and extracting the range by taking the 
square root of the fitted range -squared polynomial, the accuracy of the 
range is greatly increased as compared to fitting the range directly. 

For example, in the case of a satellite in a 1000 kilometer circular or- 
bit at the overhead portion of the orbit, for a 200 second data arc where 
the data rate is l/second and there is no random noise, bias, or peri- 
odic noise on the range data, the standard deviation of fit (between the 
square root of the fitted range-squared-polynomial and the true range 
function) is 0. 8 centimeters for a 4th degree fit to the square of range. 
For a 3rd degree polynomial fit to the square of range the standard de- 
viation of fit is 9. 09 meters. For a 4th degree polynomial fit to the 
range, the standard deviation of fit (between the fitted range polynomial 
and the true range function) is 108 meters. In order to achieve a stand- 
ard deviation of 0.8 centimeters for a least squares polynomial fit to 
range directly, a 10th degree polynomial is required. For a typical 
S-band system, range resolution is on the order of 1 to 2 meters. The 
advantage of a lower degree polynomial fit to range squared occurs be- 
cause the computation time for fitting a 4th degree polynomial to range 
squared is much less than the computation time required to fit a 10th 
degree polynomial to range directly. 

4. By multiplying the range rate by the range and fitting to this mixed data 
type, and extracting the range rate by dividing the fitted polynomial by 
the range (where the range has been previously determined by fitting a 
4th degree least squares polynomial to range squared as discussed in 
#3 on page 42), the accuracy of the range rate is greatly increased as 
compared to fitting to the range rate directly for the same degree poly- 
nomial. For example, in the case of a satellite in a 1000 kilometer cir- 
cular orbit at the overhead portion of the orbit for a 200 second data arc 
where the data rate is l/second and there is no random noise, bias, or 
periodic noise on the range rate data, the standard deviation of fit (be- 
tween the approximating function for range rate just discussed above 
and the true range rate function) is 0. 76 x 10" 3 meters /second for a 3rd 
degree least squares polynomial fit to the range times the range rate. 
For a 3rd degree polynomial fit to the range rate the standard deviation 
of fit (between the fitted range rate polynomial and the true range rate 
function) is 10. 9 meters/second. In order to achieve a standard devia- 
tion of 0. 76 x 10 _3 meters/second for a least squares polynomial fit to 
range rate directly, a 10th degree polynomial is required. For a typ- 
ical S-band system the range rate resolution is on the order of 0. 5 x 
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10~ 2 meters /second for a sample rate of one per second. The advantage 
of a lower degree polynomial fit to the range rate times the range oc- 
curs because the computation time for fitting a 3rd degree polynomial to 
this mixed data type is much less than the computation time required to 
fit a 10th degree polynomial to range rate directly . 

5. For satellites in near-Earth circular orbits ranging from 500 kilome- 
ters to 1500 kilometers and for satellites in elliptical orbits around the 
Earth (such as IMP-G), a 4th degree least squares polynomial fit to 
range squared and a 3rd degree least squares polynomial fit to the range 
rate times the range are sufficient for achieving a typical S-band range 
resolution of 1 to 2 meters and a typical S-band range rate resolution of 

0.005 meters/second. 

6. Gaussian random noise on the range and range rate data is unaffected 
by the above process of fitting and reduction of uncertainties by the 
square root of the number of points follows just as linear least squares 
polynomial smoothing. Biases are not smoothed out but remain un- 
changed in the data. Periodic noise on the data is not removed by this 
processing technique and such noise has a pronounced effect on the re- 
sults. For a fixed amplitude the greater the number of periods of peri- 
odic noise in the data span, the less the standard deviation of fit and the 
more this type of noise appears to behave like random noise. In addi- 
tion, for a fixed number of periods of this noise in the data span, the 
less the amplitude of the periodic noise the less the effect in the smooth- 
ing process and the less the standard deviation of fit. It must be em- 
phasized that the effects of periodic noise are not increased by the 
recommended process. 


10.0 RECOMMENDATIONS 

For satellites in circular Earth orbits at altitudes ranging from 500 kilo- 
meters to 1500 kilometers and for satellites in highly elliptical Earth orbits 
with low perigee, it is recommended that: 

1. A 4th degree least squares polynomial be fit to the square of the range 
measurements and the value for range obtained by taking the square 
root of the range -squared polynomial. 

2. A 3rd degree least squares polynomial be fit to the range times the 
range rate measurements (where the values of range are obtained by 
the method outlined in #1 and properly time tagged to coincide with 
the time of the range rate measurements). The value for range rate 
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is obtained by dividing the fitted polynomial by the range (where again 
the values of range are obtained by the method outlined in #1 and 
properly time tagged to coincide with the time extracted from the 
fitted rR polynomial). 

These methods will be tested in the near future by the authors on actual 
tracking data , and comparisons will be made of computation time and estimates 
of errors between the present smoothing methods and the methods proposed 
above. Estimating error in smoothing actual tracking data will be more diffi- 
cult than for this documented case of simulated data, but actual orbital uncer- 
tainties can be used for this judgement. 
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APPENDIX A 


ILLUSTRATION OF SMOOTHING ERRORS FROM FITTING 
TRACKING DATA IN THE LEAST SQUARES SENSE TO 
POLYNOMIALS OF 1st THROUGH 10th DEGREE 


A-l 
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A-2 


Figure A-l . Error (calculated polynomial minus true function) from using 1st degree polynomial as approxi- 
mation to 44 points of simulated range data at one per second data rate, a - 7378 km, e - U, 
zero time at spacecraft crossing station zenith (overhead pass). 
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Figure A-2. Error (calculated polynomial minus true function) from using 2nd degree polynomial as approxi- 
mation to 44 points of simulated range data at one per second data rate, a = 7378 km, e = 0, 
zero time at spacecraft crossing station zenith (overhead pass). 



Figure A-3. Error (calculated polynomial minus true function) from using 3rd degree polynomial as approxi- 
mation to 44 points of simulated range data at one per second data rate, a = 7378 km, e = 0, 
zero time at spacecraft crossing station zenith (overhead pass). 
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A-7 


Figure A-6. Error (calculated polynomial minus true function) from using 6th degree polynomial as approxi- 
mation to 44 points of simulated range data at one per second data rate, a = 7378 km, e = 0, 
zero time at spacecraft crossing station zenith (overhead pass) . 
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Figure A-8. Error (calculated polynomial minus true function) from using 8th degree polynomial as approxi- 
mation to 44 points of simulated range data at one per second data rate, a = 7378km, e = 0, 
zero time at spacecraft crossing station zenith (overhead pass) . 



SH3I3H MI 


A-10 


Figure A-9. Error (calculated polynomial minus true function) from using 9th degree polynomial as approxi- 
mation to 44 points of simulated range data at one per second data rate, a = 7378km, e = 0, 
zero time at spacecraft crossing station zenith (overhead pass) . 
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Figure A-10. Error (calculated polynomial minus true function) from using 10th degree polynomial as approxi- 
mation to 44 points of simulated range data at one per second data rate, a = 7378 km, e = 0, 
zero time at spacecraft crossing station zenith (overhead pass). 
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Figure A— 1 1 . Error (calculated polynomial minus true function) from using 1st degree polynomial as approxi- 
mation to 44 points of simulated range rate data at one per second data rate, a = 7378 km, 
e = 0, zero time at spacecraft crossing station zenith (overhead pass). 
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Figure A-12. Error (calculated polynomial minus true function) from using 2nd degree polynomial as approxi- 
mation to 44 points of simulated range rate data at one per second data rate, a = 7378 km, 
e = 0, zero time at spacecraft crossing station zenith (overhead pass) . 




dNOoas/sOTiaw ni hoots 


Figure A-13. Error (calculated polynomial minus true function) from using 3rd degree polynomial as approxi- 
mation to 44 points of simulated range rate data at one per second data rate, a - 73/8 km, 
e = 0, zero time at spacecraft crossing station zenith (overhead pass). 



008 



0NDD3S/Smaw NX 30333 


'x 

2 E 

D 00 


•'t 

< 


D 

DJ 


A-15 



aNODas/s^aiaw m aonaa 


A-16 


Figure A-15. Error (calculated polynomial minus true function) from using 5th degree polynomial as approxi- 
mation to 44 points of simulated range rate data at one per second data rate, a = 7378km, 
e = 0, zero time at spacecraft crossing station zenith (overhead pass). 
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igure A-16. Error (calculated polynomial minus true function) from using 6th degree polynomial as approxi 
motion to 44 points of simulated range rate data at one per second data rate, a = 7378km, 
e = 0, zero time at spacecraft crossing station zenith (overhead pass). 
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Figure A-17. Error (calculated polynomial minus true function) from using 7th degree polynomial as approxi- 
mation to 44 points of simulated range rate data at one per second data rate, a = 7378 km, 
e = 0, zero time at spacecraft crossing station zenith (overhead pass) . 
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Figure A-18. Error (calculated polynomial minus true function) from using 8th degree polynomial as approxi 
motion to 44 points of simulated range rate data at one per second data rate, a = 7378 km, 
e = 0, zero time at spacecraft crossing station zenith (overhead pass). 
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Figure A-19. Error (calculated polynomial minus true function) from using 9th degree polynomial as approxi 
motion to 44 points of simulated range rate data at one per second data rate, a = 7378 km, 
e = 0, zero time at spacecraft crossing station zenith (overhead pass). 
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Figure A-20. Error (calculated polynomial minus true function) from using 10th degree polynomial as approxi- 
mation to 44 points of simulated range rate data at one per second data rate, a = 7378 km, 
e = 0, zero time at spacecraft crossing station zenith (overhead pass). 
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Figure A-21 . Error (calculated polynomial minus true function) from using 1st degree polynomial as approxi- 
mation to 44 points of simulated elevation data at one per second data rate, a = 7378 km, 
e = 0, zero time at spacecraft crossing station zenith (overhead pass). 
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gure A-22. Error (calculated polynomial minus true function) from using 2nd degree polynomial as approxi- 
mation to 44 points of simulated elevation data at one per second data rate, a = 7378 km, 
e = 0, zero time at spacecraft crossing station zenith (overhead pass). 



saaaoaa ni aoaaa 


A -24 


Figure A-23. Error (calculated polynomial minus true function) from using 3rd degree polynomial as approxi 
motion to 44 points of simulated elevation data at one per second data rate, a = 7378 km, 
e = 0, zero time at spacecraft crossing station zenith (overhead pass). 
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Figure A-24. Error (calculated polynomial minus true function) from using 4th degree polynomial as approxi- 
mation to 44 points of simulated elevation data at one per second data rate, a = 7378 km, 
e = 0, zero time at spacecraft crossing station zenith (overhead pass) . 
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Figure A-25. Error (calculated polynomial minus true function) from using 5th degree polynomial as approxi- 
mation to 44 points of simulated elevation data at one per second data rate, a = 7378 km, 
e = 0, zero time at spacecraft crossing station zenith (overhead pass). 
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Figure A-2 7 . Error (calculated polynomial minus true function) from using 7th degree polynomial as approxi 
motion to 44 points of simulated elevation data at one per second data rate, a = 7378 km, 
e = 0, zero time at spacecraft crossing station zenith (overhead pass). 
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Figure A-28. Error (calculated polynomial minus true function) from using 8th degree polynomial as approxi- 
mation to 44 points of simulated elevation data at one per second data rate, a = 7378 km, 
e = 0, zero time at spacecraft crossing station zenith (overhead pass). 
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Figure A-29. Error (calculated polynomial minus true function) from using 9th degree polynomial as approxi 
motion to 44 points of simulated elevation data at one per second data rate, a = 7378 km, 
e = 0, zero time at spacecraft crossing station zenith (overhead pass). 
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Figure A-30. Error (calculated polynomial minus true function) from using 10th degree polynomial as approxi- 
mation to 44 points of simulated elevation data at one per second data rate, a = 7378 km, 
e = 0, zero time at spacecraft crossing station zenith (overhead pass). 
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Figure A-31 . Error (calculated polynomial minus true function) from using 1st degree, least squares fitted 
polynomial as approximation to 200 points of simulated range data at one per second data 
rate, a = 7378km, e = 0, zero time at spacecraft crossing station zenith (overhead pass). 
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Figure A-32. Error (calculated polynomial minus true function) from using 2nd degree, least squares fitted 
polynomial as approximation to 200 points of simulated range data at one per second data 
rate, a = 7378 km, e = 0, zero time at spacecraft crossing station zenith (overhead pass) . 
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Figure A-33. Error (calculated polynomial minus true function) from using 3rd degree, least squares fitted 
polynomial as approximation to 200 points of simulated range data at one per second data 
rate, a = 7378 km, e = 0, zero time at spacecraft crossing station zenith (overhead pass) . 
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Figure A-34. Error (calculated polynomial minus true function) from using 4th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range data at one per second data 
rate, a = 7378km, e = 0, zero time at spacecraft crossing station zenith (overhead pass). 
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Figure A-35. Error (calculated polynomial minus true function) from using 5th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range data at one per second data 
rate, a = 7378km, e = 0, zero time at spacecraft crossing station zenith (overhead pass). 
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Figure A-36. Error (calculated polynomial minus true function) from using 6th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range data at one per second data 
rate, a = 7378 km, e = 0, zero time at spacecraft crossing station zenith (overhead pass) . 
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Figure A-37. Error (calculated polynomial minus true function) from using 7th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range data at one per second data 
rate, a = 7378km, e = 0, zero time at spacecraft crossing station zenith (overhead pass). 
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Figure A-38. Error (calculated polynomial minus true function) from using 8th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range data at one per second data 
rate, a = 7378km, e = 0, zero time at spacecraft crossing station zenith (overhead pass). 
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Figure A-40. Error (calculated polynomial minus true function) from using 10th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range data at one per second data 
rate, a = 7378km, e = 0, zero time at spacecraft crossing station zenith (overhead pass). 
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Figure A-41 . Error (calculated polynomial minus true function) from using 1st degree, least squares fitted 
polynomial as approximation to 200 points of simulated range rate data at one per second 
data rate, a = 7378 km, e = 0, zero time at spacecraft crossing station zenith (overhead 
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Figure A-42. Error (calculated polynomial minus true function) from using 2nd degree, least squares fitted 
polynomial as approximation to 200 points of simulated range rate data at one per second 
data rate, a = 7378 km, e = 0, zero time at spacecraft crossing station zenith (overhead 
pass) . 
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Figure A-43. Error (calculated polynomial minus true function) from using 3rd degree, least squares fitted 
polynomial as approximation to 200 points of simulated range rate data at one per second 
data rate, a = 7378 km, e = 0, zero time at spacecraft crossing station zenith (overhead 
pass) . 
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Figure A-44. Error (calculated polynomial minus true function) from using 4th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range rate data at one per second 
data rate, a = 7378 km, e = 0, zero time at spacecraft crossing station zenith (overhead 
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Figure A-45. Error (calculated polynomial minus true function) from using 5th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range rate data at one per second 
data rate, a = 7378 km, e = 0, zero time at spacecraft crossing station zenith (overhead 
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Figure A-46. Error (calculated polynomial minus true function) from using 6th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range rate data at one per second 
data rate, a = 7378 km, e = 0, zero time at spacecraft crossing station zenith (overhead 
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Figure A-47. Error (calculated polynomial minus true function) from using 7th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range rate data at one per second 
data rate, a = 7378 km, e = 0, zero time at spacecraft crossing station zenith (overhead 
pass) . 




Figure A-48. Error (calculated polynomial minus true function) from using 8th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range rate data at one per second 
data rate, a = 7378 km, e = 0, zero time at spacecraft crossing station zenith (overhead 
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Figure A-49. Error (calculated polynomial minus true function) from using 9th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range rate data at one per second 
data rate, a = 7378 km, e = 0, zero time at spacecraft crossing station zenith (overhead 
pass) . 
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Figure A-50. Error (calculated polynomial minus true function) from using 10th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range rate data at one per second 
data rate, a = 7378 km, e = 0, zero time at spacecraft crossing station zenith (overhead 
pass) . 
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Figure A-51 . Error (calculated polynomial minus true function) from using 1st degree, least squares fitted 
polynomial as approximation to 200 points of simulated elevation data at one per second 
data rate, a = 7378 km, e = 0, zero time at spacecraft crossing station zenith (overhead 
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Figure A-52. Error (calculated polynomial minus true function) from using 2nd degree, least squares fitted 
polynomial as approximation to 200 points of simulated elevation data at one per second 
data rate, a = 7378km, e = 0, zero time at spacecraft crossing station zenith (overhead 
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Figure A-54. Error (calculated polynomial minus true function) from using 4th degree, least squares fitted 
polynomial as approximation to 200 points of simulated elevation data at one per second 
data rate, a = 7378 km, e = 0, zero time at spacecraft crossing station zenith (overhead 
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Figure A-55. Error (calculated polynomial minus true function) from using 5th degree, least squares fitted 
polynomial as approximation to 200 points of simulated elevation data at one per second 
data rate, a = 7378km, e = 0, zero time at spacecraft crossing station zenith (overhead 
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Figure A-56. Error (calculated polynomial minus true function) from using 6th degree, least squares fitted 
polynomial as approximation to 200 points of simulated elevation data at one per second 
data rate, a = 7378 km, e = 0, zero time at spacecraft crossing station zenith (overhead 
pass) . 
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Figure A-58. Error (calculated polynomial minus true function) from using 8th degree, least squares fitted 
polynomial as approximation to 200 points of simulated elevation data at one per second 
data rate, a = 7378 km, e = 0, zero time at spacecraft crossing station zenith (overhead 
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Figure A-59. Error (calculated polynomial minus true function) from using 9th degree, least squares fitted 
polynomial as approximation to 200 points of simulated elevation data at one per second 
data rate, a = 7378 km, e = 0, zero time at spacecraft crossing station zenith (overhead 
pass) . 
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Figure A-60. Error (calculated polynomial minus true function) from using 10th degree, least squares fitted 
polynomial as approximation to 200 points of simulated elevation data at one per second 
data rate, a = 7378 km, e = 0, zero time at Spacecraft crossing station zenith (overhead 
pass) . 
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Figure A-61 . Error (calculated polynomial minus true function) from using 1st degree, least squares fitted 
polynomial as approximation to 44 points of simulated range data at one per second data 
rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (overhead 
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Figure A-62. Error (calculated polynomial minus true function) from using 2nd degree, least squares fitted 
polynomial as approximation to 44 points of simulated range data at one per second data 
rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (overhead 
pass) . 
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Figure A-63. Error (calculated polynomial minus true function) from using 3rd degree, least squares fitted 
polynomial as approximation to 44 points of simulated range data at one per second data 
rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (overhead 
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Figure A-64. Error (calculated polynomial minus true function) from using 4th degree, least squares fitted 
polynomial as approximation to 44 points of simulated range data at one per second data 
rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (overhead 
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Figure A-65. Error (calculated polynomial minus true function) from using 5th degree, least squares fitted 
polynomial as approximation to 44 points of simulated range xlata at one per second data 
rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (overhead 
pass) . 
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Figure A-66. Error (calculated polynomial minus true function) from using 6th degree, least squares fitted 
polynomial as approximation to 44 points of simulated range data at one per second data 
rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (overhead 
pass) . 
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Figure A-6 7. Error (calculated polynomial minus true function) from using 7th degree, least squares fitted 
polynomial as approximation to 44 points of simulated range data at one per second data 
rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (overhead 
pass) . 
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Figure A-68. Error (calculated polynomial minus true function) from using 8th degree, least squares fitted 
polynomial as approximation to 44 points of simulated range data at one per second data 
rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (overhead 
pass) . 
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Figure A-69. Error (calculated polynomial minus true function) from using 9th degree, least squares fitted 
polynomial as approximation to 44 points of simulated range data at one per second data 
rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (overhead 
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Figure A-70. Error (calculated polynomial minus true function) from using 10th degree, least squares fitted 
polynomial as approximation to 44 points of simulated range data at one per second data 
rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (overhead 
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Figure A-71 . Error (calculated polynomial minus true function) from using 1st degree, least squares fitted 
polynomial as approximation to 44 points of simulated range rate data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A-72. Error (calculated polynomial minus true function) from using 2nd degree, least squares fitted 
polynomial as approximation to 44 points of simulated range rate data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
head pass) . 
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Figure A-73. Error (calculated polynomial minus true function) from using 3rd degree, least squares fitted 
polynomial as approximation to 44 points of simulated range rate data at one per second 
data rate, a = 95,000 km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A-74. Error (calculated polynomial minus true function) from using 4th degree, least squares fitted 
polynomial as approximation to 44 points of simulated range rate data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A -75. Error (calculated polynomial minus true function) from using 5th degree, least squares fitted 
polynomial as approximation to 44 points of simulated range rate data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A-79. Error (calculated polynomial minus true function) from using 9th degree, least squares fitted 
polynomial as approximation to 44 points of simulated range rate data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
head pass) . 
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Figure A-80. Error (calculated polynomial minus true function) from using 10th degree, least squares fitted 
polynomial as approximation to 44 points of simulated range rate data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A-81 . Error (calculated polynomial minus true function) from using 1st degree, least squares fitted 
polynomial as approximation to 44 points of simulated elevation data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A-82. Error (calculated polynomial minus true function) from using 2nd degree, least squares fitted 
polynomial as approximation to 44 points of simulated elevation data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A-83. Error (calculated polynomial minus true function) from using 3rd degree, least squares fitted 
polynomial as approximation to 44 points of simulated elevation data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A-84. Error (calculated polynomial minus true function) from using 4th degree, least squares fitted 
polynomial as approximation to 44 points of simulated elevation data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
head pass) . 
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Figure A-86. Error (calculated polynomial minus true function) from using 6th degree, least squares fitted 
polynomial as approximation to 44 points of simulated elevation data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A-87. Error (calculated polynomial minus true function) from using 7th degree, least squares fitted 
polynomial as approximation to 44 points of simulated elevation data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A-88. Error (calculated polynomial minus true function) from using 8th degree, least squares fitted 
polynomial as approximation to 44 points of simulated elevation data at one per second 
data rate, a = 95,000 km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A-89. Error (calculated polynomial minus true function) from using 9th degree, least squares fitted 
polynomial as approximation to 44 points of simulated elevation data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A-90. Error (calculated polynomial minus true function) from using 10th degree, least squares fitted 
polynomial as approximation to 44 points of simulated elevation data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A-91 . Error (calculated polynomial minus true function) from using 1st degree, least squares fitted 
polynomial as approximation to 200 points of simulated range data at one per second data 
rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (overhead 
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Figure A-92. Error (calculated polynomial minus true function) from using 2nd degree, least squares fitted 
polynomial as approximation to 200 points of simulated range data at one per second data 
rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (overhead 
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Figure A-93. Error (calculated polynomial minus true function) from using 3rd degree, least squares fitted 
polynomial as approximation to 200 points of simulated range data at one per second data 
rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (overhead 
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Figure A-94. Error (calculated polynomial minus true function) from using 4th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range data at one per second data 
rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (overhead 
pass) . 
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Figure A-95. Error (calculated polynomial minus true function) from using 5th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range data at one per second data 
rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (overhead 
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Figure A-96. Error (calculated polynomial minus true function) from using 6th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range data at one per second data 
rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (overhead 
pass) . 
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Figure A-97. Error (calculated polynomial minus true function) from using 7th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range data at one per second data 
rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (overhead 
pass) . 
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Figure A-98. Error (calculated polynomial minus true function) from using 8th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range data at one per second data 
rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (overhead 
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Figure A-99. Error (calculated polynomial minus true function) from using 9th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range data at one per second data 
rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (overhead 
pass) . 
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Figure A-IOO. Error (calculated polynomial minus true function) from using 10th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range data at one per second data 
rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (overhead 
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Figure A-101 . Error (calculated polynomial minus true function) from using 1st degree, least squares fitted 
polynomial as approximation to 200 points of simulated range rate data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A-102. Error (calculated polynomial minus true function) from using 2nd degree, least squares fitted 
polynomial as approximation to 200 points of simulated range rate data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A-103. Error (calculated polynomial minus true function) from using 3rd degree, least squares fitted 
polynomial as approximation to 200 points of simulated range rate data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A-104. Error (calculated polynomial minus true function) from using 4th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range rate data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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polynomial as approximation to 200 points of simulated range rate data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over 
head pass) . 
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Figure A-106. Error (calculated polynomial minus true function) from using 6th degree, least squares fitted 
polynomial as approximdtion to 200 points of simulated range rate data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
head pass) . 
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Figure A-107. Error (calculated polynomial minus true function) from using 7th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range rate data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A-108. Error (calculated polynomial minus true function) from using 8th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range rate data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A-109. Error (calculated polynomial minus true function) from using 9th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range rate data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A-110. Error (calculated polynomial minus true function) from using 10th degree, least squares fitted 
polynomial as approximation to 200 points of simulated range rate data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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polynomial as approximation to 200 points of simulated elevation data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (ovei 
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Figure A-l 12. Error (calculated polynomial minus true function) from using 2nd degree, least squares fitted 
polynomial as approximation to 200 points of simulated elevation data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A-l 13. Error (calculated polynomial minus true function) from using 3rd degree, least squares fitted 
polynomial as approximation to 200 points of simulated elevation data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
head pass) . 
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Figure A— 114. Error (calculated polynomial minus true function) from using 4th degree, least squares fitted 
polynomial as approximation to 200 points of simulated elevation data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A— 115. Error (calculated polynomial minus true function) from using 5th degree, least squares fitted 
polynomial as approximation to 200 points of simulated elevation data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A-l 16. Error (calculated polynomial minus true function) from using 6th degree, least squares fitted 
polynomial as approximation to 200 points of simulated elevation data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 



A-118 


Figure A-117. Error (calculated polynomial minus true function) from using 7th degree, least squares fitted 
polynomial as approximation to 200 points of simulated elevation data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A-l 18. Error (calculated polynomial minus true function) from using 8th degree, least squares fitted 
polynomial as approximation to 200 points of simulated elevation data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
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Figure A-120. Error (calculated polynomial minus true function) from using 10th degree, least squares fitted 
polynomial as approximation to 200 points of simulated elevation data at one per second 
data rate, a = 95,000km, e = 0.89, zero time at spacecraft crossing station zenith (over- 
head pass) . 


